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support. 
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Adaptive  Beamforming  with  Imperfect  Arrays: 
Pattern  Effects  and  Their  Partial  Correction 


1.  INTRODUCTION 


This  report  considers  adaptive  beamforming1’2’3  with  arrays  where  some  of  the 
usual  idealization  assumptions  are  not  valid.  Specifically  the  digital,  open  loop, 
implementations  (usually  called  sampled  matrix  inversion,  SMI,  or  direct  natrix 
inversion,  DMI)  of  Applebaum’s4  "Maximum  Signal  to  Noise  plus  Interference,  SNIR, 
ratio"  criterion  and  Frost’s5  "Constrained  power  minimization”  criterion  are 


Received  for  publication  15  January  1993 

1  Compton,  R.T.,  (1988),  Adaptive  Antennas ,  Prentice  Hall. 

2  Hudson,  J.E.,  (1981),  Adaptive  Array  Principles ,  Peter  Peregrinus. 

3  Monzingo,  R.  and  Miller,  T.,  (1980),  Introduction  to  Adaptive  Arrays,  John  Wiley 
&  Sons. 

4  Applebaum,  S.P.,  (1976),  Adaptive  arrays,  IEEE  Trans.  Antennas  and  Propagat., 
Sept.,  pp.  585-598. 

5  Frost,  O.L.,  (1972),  An  algorithm  for  the  linearly  constrained  adaptive 
processing,  Proc.  IEEE,  Aug.,  pp.  926-935. 
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considered.  The  received  signals  in  the  channels  are  sampled,  a  covariance  matrix 
is  formed  and  inverted,  or  transformed  by  a  similar  operation,  and  multiplied  with  a 
steering  vector  giving  an  optimal  channel  weight  vector. 

Usually  only  the  resulting  SNIR  is  considered  as  a  quality  measure  of  the 
adaptation.  However,  especially  in  radar,  it  is  also  important  to  have  a  good 
pattern  shape  with  low  sidelobes  and  a  narrow  main  beam  to  suppress  reflected 
clutter  and  give  good  direction  estimation. 

In  the  following  sections  different  types  of  errors  in  the  antenna,  which  will 
affect  the  adaptive  processes,  are  discussed  and  simulated  results  are  shown.  A 
narrowband  antenna  system  is  assumed.  Methods  for  compensating  for  some  of  the 
errors  are  also  discussed.  The  type  of  imperfections  considered  are  inequality 
between  the  analog  receiver  modules  (gain,  phase  shift,  noise  figure),  errors  in  the 
conversion  from  the  analog  IF  signal  to  complex  baseband  signal  (IQ  imbalance, 
DC-offset),  errors  in  the  conversion  from  analog  to  digital  signals  (signal 
quantization,  DC-offset),  errors  when  estimating  the  covariance  matrix  (sampling 
effects),  and  errors  in  dealing  with  the  resulting  weight  vector  (weight 
quantization).  Also,  the  effect  of  saturation  of  the  signal  and  the  corresponding 
generation  of  intermodulation  products  and  harmonics  are  considered. 

A  modification  of  Frost’s  algorithm  for  tapered  quiescent  weights  is  also 
discussed.  This  is  not  a  correction  of  an  error  but  the  modification  used  is 
similar  to  the  correction  of  some  of  the  errors  considered. 

In  the  simulations  in  the  following  sections  two  basic  configurations  have  been 
used  with  small  variations  in  some  cases.  These  basic  configurations  are: 


A) 


Linear  array  with  half  wavelength  spacing  and  /V  =  16  elements. 

O 

Desired  signal  at  broadside,  0^  =  0  . 

Uniform  quiescent  weights  (siav/.r  -  type  quiescent  pattern). 
Desired  signal  to  noise  ratio  per  channel,  SNR  =  p  Jo2  =  10  dB. 


One  jammer  close  to  the  second  sidelobe  peak,  0y  =  18  . 
Jammer  to  noise  ratio  per  channel  JNR  =  py/a2  =  40  dB. 


B)  Equal  to  configuration  A)  except  for: 

Quiescent  weights  tapered  for  40  dB  Chebyshev  pattern,  giving  an  aperture 
efficiency  of,  T|  =  0.76. 


2 


o 


Jammer  close  to  Chebyshev  pattern  sidelobe  peak,  0y  =  20 


The  optimal  gain,  SNIR  and  null  depth  without  any  errors  are: 


A)  Gain  =  11.96  dB  SNIR  =  21.96  dB  Null  depth  =  -121.4  dB 

B)  Gain  =  10.87  dB  SNIR  =  20.87  dB  Null  depth  =  -144.1  dB 


All  patterns  in  this  report  show  the  normalized  array  factor.  No  element 
patterns  are  included  in  the  plots  nor  in  the  weight  calculations. 

Figures  1  and  2  show  the  optimal  adapted  patterns  and  the  corresponding 
quiescent  patterns. 
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Basic  configuration  A 


Angle  [deg] 


Figure  1.  Adapted  -  ,  and  Quiescent  - - ,  Patterns  for  Configuration  A. 
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2.  SOME  BASIC  RELATIONS 


In  this  section  some  equations  to  be  used  in  the  following  sections  are  given. 

Let  x(r)  be  a  vector  of  dimension  N  containing  the  total  signals  in  the  N 

channels  of  the  antenna  as  a  function  of  time.  Also  let  x  ,  x.,  and  x  ,  be  the 

n  j  a 

separate  signal  vectors  of  thermal  noise,  jammers  (interference),  and  the  signal 
containing  the  desired  information  respectively,  \(t)  =  x^+  \j+  x^.  The  output 
signal  of  the  array  is:  s(t)  =  w?x(/)  where  w  is  a  vector  containing  the  channel 
weight  coefficients.  (In  this  report,  indicates  a  transpose,  a  complex 
conjugate,  and  +  a  conjugate  transpose.)  Applebaum4  has  shown  that  the  maximum 
expected  ratio  of  the  desired  signal  to  the  interference  plus  noise,  SNIR,  in  the 
output  signal  (E(|  represents  the  expectation  value  operator): 


m 


lwTx'P 


w  v  j  w  1-  wT(x.+x  )I2J 

J  y  j  rt' 


(2.1) 


is  obtained  when  the  complex  conjugate  of  the  weight  vector  is  given  by: 


w*~  $  *s  , 
a  a 


(2.2) 


where  is  the  steering  vector,  which  defines  the  desired  signal  x^  as  the  part  of 
x^  being  parallel  to  s^,  and  where  <t>A  is  the  covariance  matrix  of  the  interference 
plus  noise: 


=E((x  +  x  )(x.+  x  )+)  (2.3) 

a  j  ny  j  n 

(Often  the  covariance  matrix  is  defined  as  =  O  .) 

If  the  interference  and  noise  are  statistically  independent: 

<J>  =  E|x.x+)  +E|x  x+)  (2.4) 

a  1  j  ]'  1  n  n’ 

This  Applebaum  criterion  can  be  used  when  the  covariance  matrix  can  be  formed 
with  the  desired  signal  x^  either  absent  or  very  small,  which  is  often  true  in 
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radar.  Otherwise,  Frost’s  criterion,  in  which  (in  its  basic  form)  the  total  power  is 
minimized  subject  to  the  constraint  that  the  array  sensitivity  in  the  desired 
direction  is  constant,  can  be  used. 

m  i  n  {  I  wTx  I  J :  I  wTs  .  I  =constant  (2.5) 

w 

The  solution  for  this  is  again: 


w*~  d>F'Sj  (2.6) 

where  is  now  the  total  covariance  matrix: 

F 

<£  =  E{(x  ,+x.+x  )(x  ,+x.+x  )+)  (2.7) 

f  ’ v  a  j  n/v  d  j  n'  1  ' 

In  practice  the  expectation  value  of  the  covariance  matrixes  and  <J>p  is 
(assuming  stationary  signals)  estimated  by  a  time  average  so  that,  if  X  is  the 
signal  matrix  of  dimension  N±K  with  element  ( n,k )  equal  to  the  k01  sample  in  channel 
n,  X=[x  ,x  ,...,x  ],  we  make  the  estimation: 

12  K 

$  =  ^  XX+  (2.8) 

(The  effect  of  this  estimation  for  finite  K  is  discussed  in  Section  8.) 

Equations  (2.2)  and  (2.6)  may,  as  discussed  below,  be  manipulated  in  different 
ways,  useful  for  different  purposes. 

2.1  Desired  Signal  Absent  in  the  Covariance  Matrix 

If,  in  this  case,  there  are  no  jammers  present,  we  get  the  quiescent  Applebaum 
matrix  and  weights: 

4>  =  a2I  ;  w*~  s  .  (2.9) 

Aq  q  a 

where  a  is  the  noise  power  in  each  channel  (assumed  equal)  and  where  I  is  the 
identity  matrix.  Thus  s^  is  identical  to  the  complex  conjugate  of  the  quiescent 
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weights  of  the  array,  which  are  often  chosen  to  give  low  sidelobes.  This  is  desired 
in  radar  to  keep  clutter  and  other  undesired  off-mainbeam  reflected  signals  low. 

When  there  is  only  one  narrowband  jammer,  uncorrelated  with  the  noise,  Eq.  (2.4), 
the  matrix  inversion  formula,  [Reference  2,  Appendix  1]: 


*  - 1  +  *  - 1 

/*  +\-i  a-j  A  xx  A 

(A+xx  )  =  A 


1  +  x*A  '  'x 


(2.10) 


can  be  used  to  obtain: 


+ 

XX. 


(2.11) 


A  °  L  o2+  x*x  ,-1  °  L  a2+  p) 


J  J 


where  *  denotes  a  unit  vector  and  P  ■  is  the  total  jammer  power  P  .=  N  pj=  x*.-Xy  where 
Pj  is  the  jammer  power  per  channel.  The  weights  are  then  obtained  from: 


r . 

w  =  d)  s  ,  ~  s  ,  -  x. - — (x-s.) 

a  a  a  a  i  2  n  i  d 
J  a  +  P  ■  J 


(2.12) 


from  which  it  is  seen  that  the  pattern  is  equal  to  the  quiescent  pattern  minus  a 
jammer  cancellation  beam  with  weights  proportional  to  the  jammer  signals,  constant 
amplitude  and  linear  phase  for  a  plane  wave  jammer. 

The  gain  in  the  jammer  direction  relative  to  the  gain  in  the  signal  direction 
becomes: 

Relative  jammer  gain  = 


A+  * 
Xy'W  | 


A  +  . 

vw 


-  vv. 


vv 


A  A  A  ^  A 

X  ,X  » X  -S. 


uLl  f,.  ] 

_2  I 

n  x  .s  .  J 


Vd 


(2.13) 
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where  r ^(0 .)  =  I  xt'S^ ! 2  and  rj^(0^)  =  I  xfys^  1 2  are  the  quiescent  pattern  aperture 
efficiencies  in  the  jammer  and  desired  signal  directions. 

If  x^  =  s^,  T)  (0^)  =  1  and  the  last  term  in  the  denominator  becomes: 


A 

X 


.A  A  i.  A 

X  .x  .s  , 

d  J  J  d 


xdsd 


xcTsd 


- 


(2.14) 


If  the  jammer  is  in  the  sidelobe  region  this  last  term  is  small  so  that  in  the 

configurations  in  Section  1,  where  P.  la 2  =  16  x  104,  Eq.  (2.13)  gives  a  null  depth  of 

-104  dB  relative  to  the  quiescent  value,  T)  (0.)/ri  (0  ,),  in  agreement  with  the 

q  J  '  q  u 

simulations. 

From  Eq.  (2.13)  one  can  see  the  interesting  fact  that  the  received  jammer  power 
will  decrease  with  increasing  P.  due  to  an  improved  null  depth.  Since  the  signal 
gain,  Eq.  (2.15)  below,  will  also  decrease,  the  resulting  SNIR  will  of  course  not 
increase. 

The  resulting  gain  in  the  desired  signal  direction,  when  there  is  one  jammer 
present,  becomes: 


Signal  gain  =  N- 


|  *  +  *  I  2 

>vw  1 


IW 


=  *nq<ty 


p  ■ 
l+-i 

A  A  A  +  A 

a2 

J 

TTTrr 


i  + 


(2.15) 


where,  if  =  s^-  Eq.  (2.14)  can  be  used  to  simplify  the  equation.  For  our  basic 
configuration  A  this  gives  a  gain  decrease  of  0.08  dB,  in  agreement  with  the 
simulations. 

Sometimes  the  jammer  vector  x.  will  have  completely  random  components.  In  this 
case  the  algorithms  will  still  try  to  nullify  the  jammer  signal  but  the  resulting 


9 


•  A  * 

weight  vector  will  get  random  components,  w  -  +  Aw^.  The  mean  square  error  of 

the  weights  can  be  obtained  from  Eqs.  (2.12)  and  (A.3).  Eliminating  the  error  term 
parallel  to  s^,  which  will  not  give  any  sidelobe  contribution  (as  later  in 
Eq.  (2.25)),  gives  the  result: 


P,  ,2 


<Aw2>=  ±w\ _ _ r-  *  \--j-} 

A  ^  o2+P.(  1-1/N)  J  77(-o2+/>.J 


M  ^  a  +P .(  1-  1/N)  i 

J  J 

This  will,  for  not  too  small  N,  give  an  average  sidelobe  level  contribution  from 
Eq.  (A. 6)  of: 


(2.16) 


Sidelobe  level  = 


< Aw  > 

A 


(2.17) 


(This  assumes  of  course  that  there  is  no  error  in  itself  that  is  correlated  with 
Aw  .)  It  can  also  be  shown  that,  since  <Aw2>  is  essentially  proportional  to  the 

A  A  »  A 

projection  |XyS^|  ,  its  standard  deviation  is,  when  Xj  has  normal  distributed 
components,  approximately  equal  to  the  expectation  value  itself. 

2.2  Desired  Signal  Present  in  the  Covariance  Matrix 

This  is  the  Frost  criterion,  and  with  no  jammer  we  get,  using  the  inversion 
formula  of  Eq.  (2.10): 


*  '  .  !Ji.  .Ji?  JL 1  ■ 


n  *  *  ♦ 

P  jX  ,X_n 


(2.18) 


*  _  .  I  A  A  d  A  a  A 

W  =  O  S  ,  ~  S  ,  -  X  , - (x  ,  s  ,) 

f  Fa  a  a  2  p  a  a 

°  +  rd 


(2.19) 


More  generally,  for  uncorrelated  narrowband  signals: 
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0=0  +  x  iX  j 
f  Ada 


and,  using  Eq.  (2.10): 


O'  1  x  >x  j  O'1 

Q-'.ti-1- 

e  A  1  +  x>-  1  x . 

d  a  d 


w’  ~  0'lrs  ,-ax  ,1  =  w'-  aO'  x  , 

F  A  L  a  aJ  A  A  fl 

where  w*  is  obtained  using  the  normalized  steering  vector,  s^,  and  where 


(2.20) 


(2.21) 

(2.22) 


a  = 


P  ,x> 
d  d  a 


-  i 


d  d  a  d 


(2.23) 


is  a  scalar.  Defining: 


a  = 


a 


1-«(s  iv 


«  a  +  _  .  1  A 

P  ,x  .O  s  . 
a  a  a  a _ 

1  +  'VX'Av 


(2.24) 


where  Ax^x  is  the  part  of  the  error  vector  x^-s^  that  is  perpendicular  to  s^,  Ax^x 
=  -  sj(s^'X^).  we  can  also  write: 


WF  -  [v“'4v]  ;  WA  a'<t’*lAV  <225) 

A  A 

It  is  seen  that  w  -  w  if,  and  only  if,  x  ,  is  parallel  to  s^,  (x^  =  s^).  (The 
Frost  case  with  Xj  =  s ^  is  equivalent  with  the  Applebaum  case  with  x.  =  x^  =  m 
which  case  Eq.  (2.15)  gives  no  signal  gain  decrease.) 

Furthermore  since  a',  Eq.  (2.24),  is  often  large,  the  Frost  algorithm  is  in 
general  much  more  sensitive  to  errors  (when  x^  *  s^)  than  Applebaum’s  algorithm. 

When  Pj  is  very  large  the  second  term  in  the  denominator  of  Eqs.  (2.23)  and  (2.24) 
cannot  be  neglected  and  we  get  a  completely  error  dominated  pattern,  independent  of 
the  signal  power  Pj. 
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Also  when  there  are  errors  the  quiescent  Applebaum  aperture  efficiency, 
1^(0^)  =  |x/sj|2’  not  be  equal  to  one  and  the  Frost  gain  in  the  signal 
direction,  x^,  will  then,  according  to  Eq.  (2.15)  with  Xj  =  *  s^,  decrease  to: 


Signal  gain  =  N 


I  w 


wvv 


1  +  An 

q 


(2.26) 


where  Ar)^  =  l-ti^G^). 

As  noted  above,  the  Frost  and  Applebaum  weights,  and  hence  also  the  null  depths, 
are  identical  when  x^  =  s^.  When  this  is  not  true,  the  jammer  null  depth,  and  also 
the  null  depths  when  there  are  more  than  one  jammer,  can  be  found  by  expressing  the 
weight  vector  w  in  terms  of  the  old  weight  vector  w  .  Thus  if  w  is  the  weight 
vector,  and  the  covauance  matrix,  for  a  specific  jammer  configuration,  including 
the  desired  signal  in  the  Frost  case,  the  weights  with  one  more  jammer  present  can, 
using  Eq.  (2.10),  be  found  from: 


w*=  [i - l 

*-  1  + 


,  .  |  A  A  i 

P  <I>  x.x  . 
0  J  J 


A  a  ]  A 

P  .x  <t>  X. 

J  J  0  J 


K 


(2.27) 


which  again  gives  the  old  pattern  plus  a  jammer  cancellation  beam. 

For  one  jammer  in  the  Frost  case  is  the  covariance  matrix  given  by  Eq.  (2.18) 
and  we  get  the  jammer  gain  relative  to  the  signal  direction  gain,  equivalent  to 
Eq.  (2.13)  for  the  Applebaum  case: 

|x  .  w  | 

Relative  jammer  gain  =  — ^ — - — 


VV 

vv 


1+ 


A  ^  -A  A  A  A  A 

pd 


X  .s  , 
J  d 


P  . 

i+-4 


A  A  -A  A  ^  A  A  ,a~  A 

A  +  A 

Vrf 


(2.28) 
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which  differs  from  Eq.  (2.13)  only  in  the  numerator  containing  the  desired  signal 
power  P,.  Again  we  see  that  the  Frost  algorithm  is  quite  sensitive  to  steering 

A  A 

vector  errors,  x^  *  s^. 

For  more  jammers  the  relative  "jammer  gain  decrease"  becomes,  assuming  lwl  =  lw  I: 

I  A+  *  I  2 

lx.-W  I  J 

Jammer  gain  decrease  =  -  =  - — — —  (2.29) 

x  -w  fl  +  P  \  <&'  x.j 
jo  L  j  j  o  yJ 

which,  for  not  too  many  strong  jammers  (when  x  x.  may  be  small  even  when  Xj  is 
not  close  to  an  old  jammer),  will  still  give  a  deep  null,  keeping  the  old  ones. 

2.3  Eigenvalue  Decomposition  of  the  Covariance  Matrix 

The  properties  of  the  covariance  matrixes  and  weight  vectors  can  also  be 
analyzed  in  terms  of  eigenvalue  decomposition.  Since  the  covariance  matrixes  are 
Hermitian  there  exist  N  orthogonal  eigenvectors  e(...eN,  forming  a  unitary  matrix,  U 
=  so  that: 


UU+  =  U+U  =  I 

(2.30) 

ch  =  U  A  IT 

(2.31) 

cD'1  =  UA-'lT 

(2.32) 

where  A  is  a  diagonal  matrix  containing  the  eigenvalues,  A.  ,  of  <!>.  If  there  are  M 
uncorrelated  jammers,  including  the  desired  signal  in  Frost’s  case,  and  if  all 
channels  have  identical  noise,  there  will  be  N-M  eigenvectors  that  are  orthogonal  to 
all  signals  present,  called  noise  eigenvectors  and  spanning  the  noise  subspace,  and 
M  eigenvectors  corresponding  to  the  signals,  spanning  the  signal  subspace.6  (If  the 
noise  is  not  identical  in  each  of  the  N  channels  the  noise  covariance  matrix,  <I>  is 
not  equal  to  a  I  and  the  above  statement  is  not  strictly  true.  If  is  known. 


6  Schmidt,  R.,  (1986),  Multiple  emitter  location  and  signal  parameter  estimation, 
IEEE  Trans.  Antennas  and  Propagat.,  March,  pp,  276-280. 
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however,  the  generalized  eigenvalues  and  eigenvectors,  obtained  from  4>e=X<J>  e , 
should  be  used.  )  The  noise  eigenvalues  are  all,  or  should  all  be,  equal  to  a  . 

As  the  jammers  are  in  general  not  orthogonal,  xynxyrn,t  0»  the  eigenvectors  are 
not  identical  to  the  jammer  signal  vectors,  but  if  they  are  orthogonal  the 
eigenvalues  are  Xj=Pj+a2  and  the  jammer  vectors  are  parallel  to  the  eigenvectors. 
This  is  also  approximately  true  if  each  jammer  is  in  the  sidelobe  region  of  the 
desired  signal  and  the  other  jammers. 

Using  the  unitary  property  of  the  eigenvector  matrix  U  and  the  fact  that  all 
noise  eigenvalues  are  equal,  we  can  write: 


<r>1  =  ua‘u+  =f  rW 

u  i  i  i 

i=l 


which  gives  the  quiescent  pattern  and  signal  eigenvector  cancellation  patterns. 

Eq.  (2.35)  is  equivalent  to  the  method  used  by  Gabriel7,  Eq.  (3). 

For  one  jammer  in  Applebaum’s  and  no  jammer  in  Frost’s  algorithm,  M  =  1  and 
A.j=  P j+  o2  and  Eq.  (2.35)  is  identical  to  Eqs.  (2.12)  and  (2.19). 

When  the  desired  signal  is  present  during  adaptation  one  eigenvector,  say  e  ,  is 
approximately  parallel  to  the  steering  vector  s^,  and  term  M  in  the  sum  will  almost 
cancel  the  first  term.  In  this  case,  using: 


e  e  + 

M  M 


-  I  e.e* 

i*M  '  ‘ 


A 


I  + 


M 


+ 


7  Gabriel,  W.F.,  (1986),  Using  estimation  techniques  in  adaptive  processing 
antenna  systems,  IEEE  Trans.  Aerospace  and  Electronic  Systems ,  July,  pp.  397-401. 


(2.33) 

(2.34) 

(2.35) 


(2.36) 
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(2.37) 


M-  I  X  X  N 

wf - v.t Kvei 


i  =1  i 


0  i=M+ 1 


or,  denoting  the  error  of  s ^  being  orthogonal  to  x ^  with  As^. 
m-  i  X  X  r  m 

wf~  v  I  +[-r4HAsrf  - 

i  =  I  i  0  L  i  =  1  J 


(2.38) 


Since  X  for  i  <  M  is  usually  large,  this  is  approximately  equal  to  Eq.  (2.35) 
except  for  the  last  term. 

When  Pj,  and  hence  XM,  is  very  large  and  A is  nonzero  the  first  term,  s^,  can 
be  neglected  and  the  pattern  will  become  an  error  dominated  pattern  with  jammer 
nulls. 


2.4  Methods  of  Correction 

Some  of  the  errors  discussed  in  the  following  sections  can  be  partly  corrected 
for  by  uJng  essentially  three  different  methods.  These  are  also  discussed  more  in 
the  sections  where  they  are  used. 

In  both  the  Applebaum  and  Frost  cases  we  know  that  all  noise  eigenvalues  are,  or 
should  be  for  good  patterns,  equal  to  0  .  However,  this  is  not  always  true  for  real 
covariance  matrixes.  Therefore,  one  method  for  improvement  is  to  identify  the  noise 
eigenvalues  by  their  magnitude  and  set  them  to  this  value. 

X  -»  02  (2.39) 

i, noise 

This  is  equivalent  to  the  method  used  by  Gabriel.7 

Another  method  to  decrease  the  effects  of  errors  in  the  covariance  matrix  is 
simply  to  make  their  relative  magnitudes  smaller  by  adding  artificial  noise  to  the 
covariance  matrix,  ’  forming  an  artificial  matrix: 

cb  =  cp  +  a2  I  (2.40) 

art  art 


8  Carlson,  B.D.,  (1988),  Covariance  matrix  estimation  errors  and  diagonal  loading 
in  adaptive  arrays,  IEEE  Trans  Aerospace  and  Electronic  Systems,  July,  pp.  397-401. 
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r 


Of  course  the  noise  is  added  in  the  covariance  matrix  only  and  not  to  the  real 
signals.  Since  the  error  effects  in  the  Frost  algorithm  are  often  functions  of  P^/ a2 
this  artificial  noise  will  decrease  these  errors  too. 

In  the  Frost  algorithm  we  often  want  the  steering  vector  to  be  parallel  to 
the  desired  signal  vector  x^.  Even  if  this  vector  cannot  be  directly  found  from  the 
covariance  matrix,  it  is,  when  all  channels  have  identical  noise  (and  approximately 
when  the  noise  is  nearly  the  same),  contained  in  the  signal  subspace.  Thus, 

9  10 

modifying  the  steering  vector  by  projecting  it  onto  this  signal  subspace,:  ’ 

%  =  I  Kv*.  <2-41> 

s  i  gna  I 
$ud  space 


which  is  equivalent  to  setting  the  inverses  of  the  eigenvalues  in  Eq.  (2.33)  equal  to 
zero,  will  often  give  a  big  improvement. 


3.  NONUNIFORM  NOISE  FIGURES 


When  deriving  many  of  the  equations  in  Section  2  it  was  assumed  that  the  noise 
in  all  channels  was  equal  so  that  the  noise  part  of  the  covariance  matrix  became  the 
noise  power  times  an  identity  matrix.  If  the  noise  figures  in  the  channels  are  not 
equal  this  is  not  true  but  we  get: 


[a2  1 

fl/a>  2 

>a20 

<t>  ‘  = 

'  1/a2  0 

0  V 

NJ 

n 

0  _2 

L  1/aJ 

(3.1) 


So  that,  if  a2  =  a2  ±  A  a2  where  Aa2  is  the  deviation  from  the  mean  value,  the 

I 

quiescent  weights  will  get  the  same  relative  error. 


9  Steyskal,  H.,  (1991),  Array  error  effects  in  adaptive  beamforming.  Microwave 
Journal ,  Sept.,  pp.  101-112. 

10  Bull,  J.F.,  Amao,  M.A.,  and  Burgess,  L.R.,  (1990),  Hypersensitivity  effects  in 
adaptive  antenna  arrays,  IEEE  AP-S  Symposium ,  1990,  pp.  396-399. 
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For  rectangular  distributed  noise,  with  noise  Figures  F  ±  AF^,  the  expected 
mean  sidelobe  level  contribution  with  no  jammers  will  be,  from  Eq.  (A. 6): 


Sidelobe  level  =  -i-  -  = 

*  A 


(3.2) 


(provided  all  channels  see  290K  external  noise). 

Figure  3  shows  simulated  patterns,  for  configuration  B  using  Applebaum’s 
algorithm,  with  relative  noise  Figure  deviations  A F  =  0,  ±0.1,  ±0.2  and  ±0.4  dB 

max 

respectively.  The  expected  mean  sidelobe  contribution  level  is,  from  Eq.  (3.2), 

-48,  -42  and  -36  dB  respectively. 

If  a  is  known  this  can  be  compensated  for  to  give  good  quiescent  patterns  by 

adding  a  matrix  o2 1  -  or,  as  in  Reference  4,  use  the  modified  steering  vector 

s"J  =  <I>  s, 
a  n  d 

One  method  to  improve  the  result,  in  which  the  noise  Figures  do  not  have 
to  be  accurately  known,  is  to  make  their  relative  variation  smaller  by  artificially 
increasing  the  noise  in  the  covariance  matrix  according  to  Eq.  (2.40).  This  will 
make  the  relative  variation,  and  thus  the  sidelobe  contribution,  a  factor  1  + 
a^/a2  smaller. 

Another  method  is  to  modify  the  noise  eigenvalues  of  the  covariance  matrix.  In 
the  ideal,  uniform  noise  Figure  case,  these  are  all  equal  to  o2  so  that  identifying 
them  by  their  magnitude  and  setting  them  all  equal  to  a2,  Eq.  (2.39),  will  improve 
the  result.  As  discussed  in  Section  2  and  Reference  6  the  nonuniform  noise  Figures  will 


also  modify  the  signal  eigenvectors  slightly  but  this  effect  is  of  minor  importance. 

Figure  4  shows  patterns  when  these  two  methods  have  been  used  for  the  AF 
if).4  dB  case. 

A  A 

If  the  steering  vector  is  parallel  to  the  desired  sig  ii  ecior,  =  x^,  the 
Frost  and  Applebaum  algorithms  give  idendcal  resui  mussed  in  Section  2. 

Their  eigenvalue  spectra,  however,  are  not  equal.  The  eigenvalue  decomposition  in 
Section  2  assumed  identical  noise  Figures.  When  these  are  not  identical  the  signal 
eigenvectors  are  still  almost  the  same  but  there  will  be  a  small  noise  subspace 
component  in  the  steering  vector,  As^  in  Eq.  (2.38),  which  will  give  most  of  the 
sidelobe  level  increase.  Thus,  identifying  the  signal  subspace  by  the  magnitude  of 
the  eigenvalues  and  projecting  the  steering  vector  on  this  subspace  to  obtain  a 
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projected  steering  vector,  according  to  Eq.  (2.41),  will  improve  the  result  for 
Frost’s  algorithm. 

Figure  5  illustrates  this  for  configuration  A  by  showing  a  pattern  using  this 
method  together  with  the  original  pattern,  both  with  A =  ±1.5  dB,  and  the 
pattern  with  no  nonuniformity. 


Same  as  old  but  rand(O)  and  0.1  =+-0.4  dB 


Figure  3.  Simulated  Patterns  for  Configuration  B  Using  Applebaum’s  Algorithm.  Noise 
figure  span:  0,  ±0.1,  ±0.2  and  ±0.4  dB. 


IS 


+-0.4  dB,  -  -  +10  noise,  +100  noise,  -  mod  lambda 


Angle  [deg] 


Figure  4.  Simulated  Patterns  for  Configuration  B  Improved  by  Adding  Artificial  Noise 
and  by  Setting  the  Noise  Eigenvalues  to  a2.  Noise  figure  span  =  ±0.4  dB. 

-  Original  pattern,  -  -  -  Artificial  noise  =  IOo  , 

—  Artificial  noise  =  100o2,  -  Noise  eigenvalues  set  to  or.  (The  two  last 

patterns  are  almost  identical). 
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Relative  gain  [dB] 


noisesp=3,  -  orig,  -  steer  v  proj,  no  error 


Angle  [deg] 


Figure  5.  Simulated  Pattern  for  Configuration  A  with  Correction  by  Projecting  the  Steering 
Vector  onto  the  Signal  Subspace,  and  the  Results  for  Uniform  Noise  Figures.  Noise 

figure  span  =  ±1.5  dB.  -  Original  pattern,  -  -  -  Frost’s  algorithm  with 

steering  vector  projection,  -  No  errors.  (The  latter  two  are  coincident  in 

the  plot.) 
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4.  FROST’S  ALGORITHM  WITH  TAPERED  QUIESCENT  WEIGHTS 


To  get  a  radiation  pattern  with  low  sidelobes  for  the  Applebaum  criterion  we 
only  have  to  use  a  steering  vector,  s^,  corresponding  to  the  desired  quiescent 
pattern.  For  Frost’s  criterion  such  a  steering  vector  does  not  solve  the  problem. 
From  Eq.(2.19)  we  get  the  quiescent  weights: 


w 


Pd+°2 


(4.1) 


where  the  projection  (x^Sj)  is  of  the  order  of  one.  The  weight  vector  is  thus  a 
linear  combination  of  the  steering  vector  and  the  signal  direction  vector  x^.  If 
the  signal  is  very  weak,  P  ^  «  Gt2,  then  wf  ~  and  the  Frost  case  has  gone  smoothly 
over  into  Applebaum’s.  Figure  6  shows  the  adapted  patterns  for  configuration  B  when 
the  signal  to  noise  ratio  per  channel,  SNR,  is  equal  to  0  dB,  -12  dB  (Pyo  =  1), 

-20  dB,  and  -30  dB  together  with  the  Applebaum  pattern,  corresponding  to  SNR  = 
dB.  The  ratio  P  /(Py-a")  is  equal  to  0.94,  0.5,  0.14,  0.016,  and  0  respectively. 

It  is  seen  that  very  low  SNR  is  required  for  a  good  pattern.  When  there  is  no 
jammer  the  broadside  gain  can  be  found  from  Eq.  (2.26)  to  be  18.4,  2.3,  0.34,  and 
0.03  dB  below  the  Applebaum  gain.  To  get  good  patterns  even  for  stronger  signals 
the  following  modifications  of  the  Frost  algorithm  may  be  used. 

A)  Increase  the  noise  artificially  by  adding  a  matrix  a  I,  Eq.  (2.40),  to  the 

covariance  matrix  so  that  P  ,  «  (f  =  d2  +C?2  .  This  will  result  in  a  slightly  lower 

a  to(  art 

SNIR  when  a  jammer  is  present  although,  of  course,  the  added  noise  is  not  present 
in  the  real  signal.  From  Eq.  (4.1)  it  is  seen  that,  with  no  jammer  present,  this  is 
equivalent  to  a  low  power  signal  as  in  Figure  6. 


B)  The  problem  with  Frost’s  criterion  is  that  the  desired  signal,  x ^  destroys  some 
nice  properties  of  the  covariance  matrix.  If  this  signal  could  be  subtracted  from 
the  total  signal  the  problem  would  be  solved.  However,  as  noted  from  Figure  6,  for 
the  Frost  case  to  become  an  Applebaum  case,  it  is  necessary  that  the  residual 
signal  power  A P ^  «  C2,  or  that  the  residual  signal-to-noise  ratio  per  channel 
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Angle  [deg] 


Figure  6.  Adapted  Patterns  Using  Frost’s  Algorithm,  Configuration  B.  Signal  to 

noise  ratio  per  channel,  SNR,  =  -  0  dB,  -  -  -  -12  dB, .  -20  dB, 

. 30  dB  and  -  -«>  dB. 


(SNR)  es  «  UN,  which  is  not  easy  to  obtain.  One  possible  method  could  be  to 
use  signals  from  subarrays  with  a  null  in  the  desired  signal  direction  to  obtain 
the  adaptive  weights  but  apply  them  to  subarrays  with  a  maximum  in  this 
direction  to  get  the  antenna  signal,  as  in  the  Duvall  beamformer.11  This  assumes  in 
principle  that  all  antenna  element  patterns  are  identical,  and  that  the  signal  and 
jammers  are  uncorrelated  plane  waves  so  that  the  optimal  pattern  has  deep 
directional  nulls. 

C)  If  the  artificial  noise  increase  is  different  for  the  channels,  so  that  the 
total  artificial  covariance  matrix  is: 

4>  =  <D  +  o2  D  -  a2I  (4.2) 

art  art 

where  D  is  a  diagonal  matrix,  the  total  noise  part  of  the  covariance  matrix 
becomes:  (  =  a^D.  Using  Eq.  (2.25)  with  d>A=  the  quiescent  weights  can 

be  written  as: 


D1 


Sd' 


Pd 


0  _  A  4.  _  .  I  A 

a  +  P  ,\  ,  D  Ax  ,, 

art  a  u  a-L 


-Av 


-  D  -  cTAx^ 


(4.3) 


where  a"  is  a  scalar. 

Thus,  if  Sj  is  the  usual  plane  wave  steering  vector,  equal  to  x^  and  if  D 
contains  the  inverse  of  the  desired  amplitude  weights  such  as  Dolph-Chebyshev,  D..  = 
1/wd.,  the  desired  quiescent  weights  are  obtained.  (If  the  artificial  noise  is 
allowed  to  be  correlated  between  the  channels,  quiescent  weights  with  irregular 
phase  may  also  be  obtained  while  keeping  <t>  Hermitian.) 

w*  ~  D"‘s  ,  (4.4) 

Fq  a 

"Widrow,  B.,  Duvall,  K.M.,  Gooch,  R.P.,  and  Newman,  W.C.,  (1982),  Signal 
cancellation  phenomena  in  adaptive  antennas:  causes  and  cures,  IEEE  Trans. 

Antennas  and  Propagat.,  May,  pp.  469-478. 
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and  in  general: 


(4.5) 


This  method  corresponds  physically  to  a  situation  where  the  antenna  edge 
elements  are,  for  tapered  quiescent  weights,  more  noisy  than  the  central  ones 
(compare  Section  3).  To  get  maximum  SNIR  the  antenna  will  then  put  more  weights  on 

A 

the  central  elements.  It  can  also  be  used  for  Applebaum’s  algorithm.  As  is  now 
parallel  to  x^  the  result  will  then,  according  to  Eq.  (2.25),  be  identical  to  that  of 
the  Frost  algorithm. 

The  relative  jammer  gain,  or  null  depth,  can  be  found  in  the  same  way  as  when  D 
=  I,  Eq.  (2.28)  with  the  result: 


Relative  jammer  gain  = 


VV 
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1+ 


P j  x*D'  (s^s^-x^x*)D" 


i+ 


_  A  4  __  .  1  -A  A  x  A  A  x .  .  I  A 

Pj  x^D  (sjS  j-x^.xj]  D 


a  r  i 


X<*D  sd 


(4.6) 


where  T)^(0y)  =  |xyD'‘s^|2/<w2>  and  1)^(0^)  =  |x^D1s^|2/<w2>  are  the  quiescent  aperture 
efficiencies  in  the  jammer  and  desired  signal  directions  and  <w2>  is  the  mean  square 
of  the  elements  in  D1. 

When  x^  =  the  equation  simplifies  to: 
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Relative  jammer  gain  (x^  =  s^) 
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<w  > 
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1  1  (0J  >) 


(4.7) 
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where  <wd>  is  the  mean  value  of  the  elements  in  D'1  and  where  now 

=  <wd>V<wd>- 

Figure  7  shows  radiation  patterns  for  basic  configuration  B  using  both  the 
Applebaum  algorithm  with  a  Chebyshev  type  steering  vector  and  the  above  method  using 
the  diagonal  elements  D. .  =  l/wDj,  where  wDj  are  the  Chebyshev  weights  with 
w  =1,  and  with  a  =  a.  The  resulting  null  depths  are  -139  dB  using  this 

Di,max  an 

method,  which  is  also  obtained  from  Eq.  (4.7),  and  -144  dB  for  the  Applebaum  pattern 
using  the  Chebyshev  steering  vector. 

When  there  are  errors  in  the  array  that  are  corrected  for  by  some  method  such  as 
a  steering  vector  projection,  Eq.  41),  or  by  noise  eigenvalue  modification,  Eq. 

(2.39),  the  covariance  matrix  mou.  ication  in  Eq.  (4.2)  should  be  made  as  the  last 
step. 
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Relative  gain  [dB] 


Frost,  and  Applebaum,  with  Chebyscheff  quiescent  pattern 


Angle  [deg] 


Figure  7.  Radiation  Patterns  for  Configuration  B.  - Modified  noise 

covariance  matrix, - Applebaum’s  Algorithm  with  Chebyshev  Steering  Vector. 
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5.  GAIN  AND  PHASE  SHIFT  ERRORS 


These  errors  can  in  principle  be  measured  by  a  calibration  routine  and  corrected 
for,  but  in  this  Section  the  case  when  this  correction  has  not  been  made  is 
considered.  They  were  also  discussed  in  Steyskal.9,12 

When  the  channels  have  errors  in  gain  and  phase  shift  the  signal  vector  into  the 
beamformer  will  be  deformed.  Since  the  adaptive  algorithm  does  not  assume  anything 
about  the  amplitude  and  phase  relations  for  the  jammers,  for  example  a  plane  wave 
assumption,  the  jammers  will  still  be  eliminated  with  essentially  the  same  factor  as 
without  these  errors.  Due  to  the  gain  and  phase  shift  errors  however,  the  steering 
vector  Sj  is  no  longer  optimum  and  an  optimal  pattern  will  not  result. 

5.1  Random  Error  Effects  in  the  Applebaum  Algorithm 

If  the  gain  and  phase  shift  errors  are  in  front  of  the  thermal  noise  sources, 
which  we  will  assume  are  all  equal,  the  noise  covariance  matrix  will  still  be  a  I 
and  the  quiescent  Applebaum  weights  will  be  identical  to  the  error  free  weights, 

Wqo~  sd.  For  random  errors  this  will  give  an  increase  in  the  average  sidelobe 
levels  and,  as  seen  from  Eqs.  (2.12)  and  (2.35),  this  increase  is,  for  the  adapted 
pattern,  essentially  the  same  as  for  the  quiescent  pattern,  described  by  an  added 
sidelobe  level  contribution,  Eq.  (A.6): 

2  2 

Sidelobe  level  =  (5.1) 

where  e  and  8  are  the  gain  and  phase  shift  RMS  errors  relative  to  their  mean  values 
respectively. 

If  the  errors  are  after  the  thermal  noise  sources  there  will  also  be  variations 
in  the  output  noise  power  in  the  channels,  as  discussed  in  Section  3,  which  is 
correlated  with  the  gain  errors.  The  covariance  matrix  is  now  given  by  <t>  =  Cd>oC+ 
where  C  is  a  diagonal  matrix  containing  the  gain  and  phase  shift  errors, 


12  Steyskal,  H.,  (1990),  Array  error  effects  in  adaptive  beamforming,  JINA  90,  Conf. 
Proc.,  pp.  475-478. 
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C.  =  exp(£+/8),  and  where  <t>Q  is  the  error  free  covariance  matrix.  By  correcting 
only  the  steering  vector,  by  setting  =  Cs^  (so  that  s^  =  x^  for  non-tapered 
steering  vectors),  the  weight  vector  becomes  w  =  C’'w  ,  which  is  the  error  free 
weights  wq  corrected  for  the  gain  and  phase  shift  errors  and  giving  the  ideal 
pattern.  When  this  steering  vector  correction  is  not  made  there  will  be  weight 
errors: 

Aw  =  <f>  \sd  -  scd) 

If  the  errors  are  small,  the  relative  errors  in  the  complex  quiescent  weights, 

A  =  Ah'  /w  ,  when  the  gain  errors  are  in  front  of  or  after  the  noise  sources  are 
v  v  , 

related  by:  A y  =  -A^,  so  that  the  added  sidelobe  level  contribution  is  still  given  by 
Eq.(5.1).  Normally  there  will  of  course  be  a  combination  of  the  two  situations 
discussed. 

Figures  8  and  9  illustrate  results  using  Applebaum’s  algorithm,  Eq.  (2.2),  with 

O 

normally  distributed  errors  with  standard  deviations  0.3  dB  and  2  ,  for 

configurations  A  and  B  respectively,  together  with  the  quiescent  patterns  and  the 

adaptive  patterns  without  any  errors.  The  sidelobe  contributions  from  Eq.  (5.1)  are 

-38  and  -37  dB  respectively.  The  null  depths  were  -122  dB  in  the  simulation  shown 

in  Figure  8,  and  -141  and  -147  dB  in  Figure  9  when  the  errors  were  in  front  of  and 

after  the  noise  sources  respectively.  The  calculated  results,  using  Eq.  (2.13)  and 

the  error  free  values  of  t)  (9  )  and  Tl  (0  ,),  were  -121  dB  and  -144  dB  respectively. 

q  J  q  0 

For  configuration  B  in  Figure  9,  where  the  induced  sidelobe  level  is  higher  than  the 
regular  sidelobes,  the  true  T) ^(€K)  may  vary  significantly  between  different 
simulations. 

Note  also  that  when  there  are  gain  errors,  the  true  antenna  aperture  efficiency 
is  not  given  by  T|  =  lx  ,-s^l2/lw  I2,  which  is  the  efficiency  referred  to  the  channel 
outputs,  but  1)  =  I  x^-s^  1 2/ 1  w+CC+w  1 2.  This  will  give  some  errors  when  using 
Eq.  (2.15)  for  the  signal  gain. 

5.2  Random  Error  Effects  in  the  Frost  Algorithm 

In  the  Frost  algorithm  the  deformed  signal  x^  is  also  included  in  the  covariance 
matrix.  If  x^  =  Sj>  as  when  correcting  the  steering  vector  as  described  above,  the 


(5.2) 
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Relative  gain  [dB] 


last  term  in  Eq.  (2.25)  will  vanish,  as  will  the  last  terms  in  Eqs.  (2.37)  and  (2.38) 
and  the  Frost  weights  will  be  identical  to  the  Applebaum  weights.  When  there  are 
uncorrected  gain  and  phase  shift  errors,  will  have  components  in  the  noise 


-  before,  —  after,  quiescent,  ..  error  free 


Angle  [dB] 


Figure  8.  Configuration  A,  Applebaum’s  Algorithm  with  0.3  dB  and  2  RMS  Errors. 

-  errors  in  front  of  the  noise  sources, - errors  after  the  noise 

sources,  - . -  quiescent  pattern,  .  error  free  pattern 
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-  error  bef,  --  errro  aft,  quiescent,  ...  errorfree 


Figure  9.  Configuration  B.  Same  parameters  as  in  Figure  8. 


0 


3' 


subspace  which,  as  X Jo2  may  be  large,  can  have  a  substantial  effect.  Comparing 
Eq.  (2.38)  with  Eq.  (2.35)  we  see  that  the  total  error  will  be  approximately  X Jo2 
larger  than  for  the  Applebaum  algorithm.  Since  XM  is  approximately  equal  to  P ,  + 
o2-  Npj  +  a ~  this  factor  can  be  very  large  for  strong  signals  and/or  many  elements. 
Alternatively,  as  a'  in  Eq.  (2.24)  is,  for  not  too  strong  signal  power,  approximately 
equal  to  P Jo2,  we  get  from  Eq.  (2.25): 


Aw 

F 


Aw  + 

A 


Aw, 


(5.4) 


where  Awa  is  the  same  error  as  in  the  Applebaum  algorithm  and  where  Aw^  is  equal  to 
Aw^  when  the  gain  errors  are  after  the  noise  sources,  remaining  essentially 
independent  of  the  position  of  these  errors  and  being  equal  to  zero  when  the 
corrected  steering  vector,  =  Cs^,  is  used.  Adding  the  two  terms  quadratically 
the  sidelobe  level  contribution  is  from  Eq.  (A. 6),  assuming  small  errors: 


Sidelobe  level  = 


(5.5) 


Physically  this  large  effect  on  the  Frost  algorithm  can  be  thought  of  as  an 
attempt  to  cancel  the  orthogonal  signal,  Xj  -  which  is  considered  as  a 

coherent  jammer  by  the  algorithm. 

Figure  10  shows  the  same  type  of  patterns  as  Figure  8  when  using  Frost’s 
algorithm,  Eq.  (2.6),  giving  a  much  worse  result  (in  this  example,  P Jo  +  1  =  161  and 
Xja  =  158).  Since  the  second  error  term  in  Eq.  (5.4)  will  now  dominate,  the 
patterns  with  erors  in  front  of  and  after  the  noise  sources  are  almost  identical. 

The  resulting  signal  gain  in  this  simulation  was  only  -4.2  dB.  With  0.3  dB  and 

O 

2  RMS  errors,  the  expected  quiescent  Applebaum  efficiency  is,  using  Eq.  (A. 2), 

T)  (9  ,)  =  0.998  and  the  Frost  gain,  using  Eq.  (2.26)  is  -5.8  dB. 
q  a 
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-  after,  —  in  front,  quiescent, ...  error  free 


Figure  10.  Configuration  A,  Frost’s  Algorithm.  Same  parameters  as  in  Figure  8. 


5.3  Improvement  of  the  Frost  Algorithm 

The  poor  result  for  Frost’s  algorithm  in  the  previous  section  can  be  improved 
however.  As  discussed  above,  the  major  cause  of  the  problem  is  that  the  signal 
vector  is  not  parallel  to  the  steering  vector,  that  is,  x^  Sj.  We  know  however 
that  the  signal  vector,  x^,  is  contained  in  the  signal  subspace  but  has  no  component 
in  the  noise  subspace  when  the  channels  have  identical  noise.  One  method  to  partly 
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cure  the  problem  is  then  to  identify  the  signal  subspace  by  the  magnitude  of  the 
eigenvalues  and  project  the  steering  vector  on  this  subspace  to  obtain  a  projected 
steering  vector,  Eq.  (2.41),  as  was  also  done  in  Section  3.  This  will  make  the  last 
term  in  Eqs.  (2.37)  and  (2.38)  equal  to  zero. 

In  the  quiescent  case  (desired  signal  but  no  jammer  present),  with  identical 
channel  noise  after  the  gain  errors,  the  projected  steering  vector  and  the  resulting 
weight  vector  will  be  exactly  parallel  to  the  signal  vector:w^~  ~  x^  =  Cx^Q, 
where  x^q  is  the  error  free  signal  vector.  Since  the  ideal  weight  vector  is,  from 
Section  5.1,  w  ~  C’x*. ,  we  get: 

ideal  <20  ° 


w  ~  C+Cw 
Fq 


ideal 


(5.6) 


2  2 

and  the  weights  will  get  a  normalized  mean  square  error  of  <A  >  =  4e  independent  of 
the  phase  error  6.  This  will  give  a  sidelobe  contribution  of: 


Sidelobe  level  =  — (5.7) 

o 

which  in  our  example  with  0.3  dB  and  2  errors  will  give  a  sidelobe  contribution 
3  dB  higher  than  for  the  Applebaum  case,  Eq.  (5.1). 

When  the  gain  errors  are  after  the  noise  sources  the  projected  steering  vector 
is  almost  identical  to  that  above,  giving  the  same  sidelobe  level.  In  Section  5.1 
we  found  that  using  a  corrected  steering  vector  =  Cs^  will  result  in  the  ideal 
pattern.  Here  we  have  indeed  =  s^.  However  since  they  are  not  exactly  equal  and 
since  the  noise  eigenvalues  are  much  smaller  than  the  signal  eigenvalues  (provided 
the  signal  is  not  very  weak)  the  last  term  in  Eq.  (2.38)  is  still  significant  when 


This  steering  vector  projection  method  assumes  of  course  that  the  signal  and 
noise  subspaces  can  be  correctly  identified.  When  this  may  be  in  doubt,  in  practice 
when  P  Jo2<  1,  the  elimination  of  the  noise  eigenvectors  is  not  so  important. 

Another  method  to  improve  the  pattern  is  to  make  the  factors  P Jo  and 
smaller  by  artificially  increasing  the  noise  in  the  covariance  matrix  to  give  a  new 
artificial  covariance  matrix  as  in  Eq.  (2.40).  The  sidelobe  contribution  is  still 
given  by  Eq.  (5.5)  with  a2  replaced  by  =  o2  +  o2^.  To  get  a  good  pattern  the 
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artificial  noise  required  is  dependent  on  the  total  power,  P^,  received  by  the 
antenna. 

Figure  11  shows  patterns  when  all  noise  subspace  components  of  the  steering 
vector  have  been  eliminated,  giving  sidelobe  contributions  of  -35  dB  from  Eq.  (5.7), 
and  patterns  with  added  noise  a  =  10  a  and  100  c r,  giving  sidelobe  contributions 

art 

of  -14  and  -31  dB  from  Eq.  (5.5).  The  signal  gain  was  12,  10,  and  12  dB  respectively 
in  agreement  with  the  gain  obtained  from  Eq.  (A.2)  with  <A2>  =  4e\  and  from 
Eq.  (2.26)  replacing  o2  with  a2o[  using  T|^(0^)  =  0.998  as  in  Section  5.2. 

When  a  low  sidelobe  pattern  is  desired,  the  method  described  in  Section  4, 

Eq.  (4.2),  can  be  used  if  the  steering  vector  projection,  Eq.  (2.41),  is  first  used  to 
correct  for  the  errors  or  if  Eq.  (4.2)  is  used  with  a  large  a2  .  For  the  steering 

art 

vector  projection  method  Eq.  (5.6)  is  still  valid,  giving  a  sidelobe  contribution  of: 

4f2 

Sidelobe  level  =  TjTr  (5'8) 

If  Eq.  (4.2)  is  used  directly  we  can  use  Eq.  (2.25)  with  =  o2  D.  Performing 

A  art 

the  same  type  of  analysis  as  for  the  Applebaum  case  we  obtain  the  sidelobe  level 
contribution,  equivalent  to  Eq.  (5.5): 


2  2  r 

Sidelobe  level  =  E-*jy  ^1  +  — 7"<wd>2]  (5.9) 

^  a  rt 

2  1 
where  <wq>  is  the  mean  value  of  the  weight  coefficients  in  D  . 

Figure  12  shows  similar  patterns  as  Figure  1 1  but  for  configuration  B.  The 

expected  sidelobe  level  contributions  from  Eqs.  (5.8)  and  (5.9)  are  -34,  -18,  and 

-34  dB  respectively. 

5.4  Null  Depth  for  the  Frost  Algorithm 

The  null  depth  for  the  Frost  algorithm  may  be  found  using  Eqs.(2.28)  or  (4.6). 

The  unit  signal  vector  can,  to  second  order  in  e  and  5,  be  written  as,  [compare 
Eq.  (8.26)]: 
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(5.10) 


where  is  the  error  free  unit  signal  vector  and  r  is  a  random  vector  where  each 
element,  except  along  x ^  ,  has  expectation  value  =  0  and  standard  deviation  =  1. 
Using  this  in  Eq.(4.6)  with  x^Q=  and  assuming  that  the  jammer  is  in  the  sidelobe 
region  we  obtain: 


Null  depth  = 


2  £2  +  52 


VV  +  ^ 

a  r  t 

vv(‘ +  -rr^<:'v) 


(5.11) 


where  t|  (0^)  and  t]  (0 .)  are  the  quiescent  (no  jammer  nor  desired 
efficiencies  and  where  r.  is  the  projection  r.  =  r+D'lx..  When  no 
projection  is  used,  I  r.  I  ^  has  the  expectation  value: 


signal)  aperture 
steering  vector 


(5.12) 


When  the  steering  vector  method  is  used  the  null  depth  is  essentially  the  same  as 
above  because  the  projection  of  the  steering  vector  on  the  jammer  vector  x .  is  not 
changed  apart  from  a  scaling  constant  close  to  one  (unless  the  errors  are  very 
large).  When  using  the  modification  of  the  covariance  matrix  in  Eq.  (4.6)  we  get, 
since  r  is  now  essentially  parallel  to  Xj,  with  expected  length  =  1: 


I  r/=  <wD>2  (5.13) 

For  configuration  A,  where  <w2>  =  <wq>  =  1,  Eq.  (5.11)  gives  -98  dB  for  the  direct 
inversion,  -98  dB  for  the  steering  vector  projection,  -97  dB  for  a2fl=  10  a2  and 
-81  dB  for  a2  =  100  a2,  while  the  results  for  configuration  B  are  -97,  -98,  -97, 

art 

and  -93  dB  respectively.  These  results  are  also  close  to  the  average  results  from 
several  simulations.  Since  the  null  depths  depend  on  one  single  projection  of  a 
random  vector  on  the  jammer  vector,  the  results  vary  significantly  between  different 
simulations. 
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5.5  Effects  of  Pointing  Error 


A  special  type  of  phase  error  occurs  when  there  is  an  antenna  pointing  error 
giving  a  linear  phase  error  in  s^,.10  This  will  have  no  effect  on  the  Applebaum 
pattern  shape  but  will  of  course  give  a  lower  gain  in  the  desired  signal  direction. 
For  Frost’s  algorithm  there  will,  as  for  the  random  errors,  be  a  stronger  effect, 
which  can  be  partly  corrected  for  in  the  same  way  as  the  random  errors. 

Figure  13  shows  patterns,  for  configuration  A  and  Frost’s  algorithm,  when  there 

O 

is  a  1.6  pointing  error  (approximately  1/4  of  the  3  dB  beam  width)  using  no 
corrections,  with  o '  =  100  c  ,  and  when  using  the  signal  subspace  steering  vector. 

The  resulting  signal  direction  gains  were  -24.6,  8.6,  and  12.0  dB  respectively.  The 
gain  values  for  the  direct  inversion  and  for  c  =  100  o  are  also  obtained  using 

art 

Eq.  (2.26).  As  seen,  the  pattern  for  the  projected  steering  vector  will  try  to 
follow  the  true  signal.  For  the  Applebaum  case,  where  the  beam  will  point  in  the 
steering  vector  direction,  the  gain  was  11.2  dB. 


6.  ERRORS  IN  QUADRATURE  SIGNAL  GENERATION 


In  analog  quadrature  signal  (IQ)  generation  the  IF-signal,  aocos((o>+(0LO)r+\j/),  is 

mixed  with  the  LO-signals  cos(u)LQ0  and  sin(a)LO0  to  generate  the  I  and  Q  video 

signals  a  =  aocos(o)r+iy)  and  a^=  aosin(GJH-t|/),  which  together  represent  the  complex 

video  signal  a  =  a  exp(j(u>/+vj/)).  There  may  however  be  an  unbalance  in  amplitude  (a) 
* 

and  phase  (P)  which,  like  the  gain  and  phase  errors  in  Section  5,  may  be  corrected 
for  if  they  are  known.14  Without  such  correction  the  signals  become: 


13  Barton,  P.,  (1980),  Digital  beam  forming  for  radar,  IEEE  Proc.,  Pt.  F,  Aug.,  pp. 
266-277. 

14  Churchill,  F.E.,  Ogar,  G.W.,  and  Thompson,  B.J.,  (1981),  The  correction  of  I  and 
Q  errors  in  a  coherent  processor,  IEEE  Trans,  on  Aerospace  and  Electronic  Systems , 
Jan.,  pp.  131-137. 
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Figure  13.  Patterns  for  Configuration  Using  Frost’s  Algorithm.  1.6  pointing 
error.  -  no  correction, - o'  =  100a2,  -  signal  subspace  steering  vector 


(6.1) 


a  i  =aQ[co  sH**  2  co  s'P- |s  i  n^] 
a^=aQ[s  in'F  -  ^  s  i  n'F-^cos1^] 

where  *¥  =  (0/+l|/.  The  complex  signal  a  can  now  be  written  as: 

ac  =  a}  +  jaQ  =  a+  a-2 $  a  =  a  +  2^  (6.2) 

where  a  is  the  correct  signal  and  where  A  is  the  total  complex  unbalance. 

If  the  video  frequency  w  *  0  the  signal  consists  of  the  correct  signal  at  w  with 
phase  \| /,  aoexp[j(cof+V}r)],  plus  an  error  signal  at  -0)  with  phase  -\j/,  2a0exPt-j(w+V)]. 
that  is,  corresponding  to  the  image  frequency  of  the  I  and  Q  mixers. 

Often  the  errors  of  the  IQ  converters  are  uncorrelated  between  the  channels. 
Furthermore,  if  the  antenna  element  modules  have  large  gain  and  phase  shift  errors, 
which  are  compensated  for  by  using  a  calibration  routine  (separating  out  the 
(D-component),  these  corrections  will  tend  to  decorrelate  otherwise  correlated  IQ 
error  signals. 

For  digital  IQ  generation,  using  a  Hilbert  filter  to  get  the  Q  signal,  the 
result  will  be  the  same  as  in  Eq.  (6.1)  with  the  difference  that  the  phase  error  p=0 
and  that  the  amplitude  error  is  given  by  the  length  of  the  filter,  and  is  not  a 
random  variable. 

6.1  Error  Effects  on  the  Quiescent  Patterns 

For  a  non-adaptive  antenna,  an  uncorrelated  error  signal  at  the  image  frequency 
will  correspond  to  an  average  sidelobe  level  contribution  in  the  total  power  pattern 
of: 


Sidelobe  level  =  |  ^  ^  (6.3) 

while,  if  the  errors  are  really  identical  in  all  channels,  the  error  signals  will 
add  to  form  a  sidelobe  peak  of: 
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Sidelobe  peak  level  =  l-  I A  1 2 


(6.4) 


Figure  14  shows  quiescent  patterns  for  configuration  B  with  Gaussian  unbalances 
having  standard  deviations  of  a=0.07  and  P=0.07,  being  equivalent  to  the  0.3  dB  and 

O 

2  amplitude  and  phase  errors  in  Section  5,  and  giving  a  sidelobe  contribution  of 
-37  dB  from  Eq.  (6.3).  The  figure  shows  the  separate  patterns  for  the  true  and  the 
image  frequencies  and  the  total  pattern  valid  for  to=0.  From  this  figure  we  can  see 
that  the  quiescent  pattern  for  the  true  frequency  is  not  affected  by  the  IQ  errors. 

6.2  Pattern  Effects  for  Nonzero  Video  Frequency 


In  the  adaptive  process  both  frequency  components,  to  and  -to,  are  included  in  the 
covariance  matrix  and  the  antenna  will  try  to  cancel  both  signals  so  that  there  will 
still  be  a  good  null  in  the  jammer  direction. 

The  image  frequency  component  of  the  jammer  will  give  a  random  signal  in  each 
channel  with  total  expected  power  P  ■< A2>/4  =  N pj< A2>/4.  Using  this  in  Eq.  (2.17)  we 
get  the  average  sidelobe  level  contribution  in  the  true  frequency  pattern: 


Sidelobe  level  = 


1 


_4oL] 

<A  2>pJ 
J 


7 


(6.5) 


Figure  15  shows  the  adapted  patterns  for  Applebaum’s  algorithm,  configuration  B, 
when  to  *  0  and  with  the  same  unbalances  as  in  Figure  14,  giving  -23  dB  sidelobes  from 
Eq.  (6.5).  This  average  sidelobe  contribution  has  a  standard  deviation  equal  to  the 
average  itself  so  that  the  average  sidelobe  level  can  vary  considerably  between 
different  simulations.  The  image  frequency  pattern  also  has  a  good  null  while  its 
average  level  relative  to  the  true  frequency  peak  is  the  same  as  the  quiescent 
pattern  sidelobe  level,  -37  dB  from  Eq.  (6.3). 

Since  the  signals  at  the  true  and  image  frequencies  are  uncorrelated  with 
totally  different  phase  relations  between  the  channels,  two  degrees  of  freedom  of 

O 

the  array  are  needed  to  get  a  null  at  20  for  both  frequencies. 

Figure  16  shows  the  same  type  of  patterns  for  configuration  A  for  both  the 
Applebaum  and  Frost  algorithms,  which,  as  seen,  give  very  similar  results.  The 
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Figure  14.  Quiescent  Patterns  for  Configuration  B.  Uncorrelated  errors 

a  =  P  =  0.07  (RMS).  -  Total  pattern  (o>=0),  - - tiue  frequency, 

- -  image  frequency. 


Relative  gain  [dB] 


IQ-err,  w  =0,  -  true  f,  -  -  image  f,  no  error:  seed=1000 


Figure  15.  Adapted  Patterns  for  Configuration  B,  Using  Applebaum’s  Algorithm. 

Uncorrelated  errors  a  =  P  =  0.07  (RMS).  -  true  frequency,  - - image 

frequency,  -  no  IQ  errors. 


Relative  gain  [dB] 


Figure  16.  Adapted  Patterns  for  Configuration  A,  (0  0,  True  and  Image  Frequencies. 

Uncorrelated  errors  a  =  P  =  0.07  (RMS).  -  Applebaum’s  algorithm, 

- Frost’s  algorithm,  -  no  IQ  errors. 


reason  is  that  when  to  *  0,  the  true  and  image  frequency  signals  are  uncorrelated  so 

that  the  error  of  the  desired  signal,  A/2  a^,  is  simply  considered  as  one  more 

uncorrelated  jammer  by  the  Frost  algorithm.  The  sidelobe  level  contribution  due  to 

this  signal  is  also  given  by  Eq.  (6.5)  with  p ■  replaced  by  p  ,,  giving  -34  dB  in  this 

J  ^ 

configuration.  One  more  degree  of  freedom  will  however  be  required. 

6.3  Pattern  Effects  with  Zero  Video  Frequency 


When  the  video  frequency  for  one  of  the  signals  is  zero  the  corresponding  error 
signal  A/2  a  is  completely  correlated  with  the  true  signal  ;  ‘here  will  not  be 
any  new  image  signal  frequency  generated  but  the  signal  vector  will  be  deformed  as 
in  Section  5. 

For  the  Applebaum  case  with  jammer  video  frequency  oiy  =  0,  but  desired  signal 
video  frequency  (0^  ^  0,  the  total  jammer  vector  is  essentially  a  plane  wave  and  the 
effect  on  the  desired  signal  pattern  is  very  small.  The  expected  sidelobe  level 
contribution  is,  as  seen  using  Eq.  (2.12),  approximately  given  by  Eq.  (6.3)  multiplied 
by  the  total  quiescent  pattern  gain  in  the  jammer  direction,  which  is  also  given  by 
Eq.  (6.3)  and  by  the  regular  sidelobe  pattern.  Close  to  the  jammer  direction 
however,  the  true  and  image  frequency  patterns  will  have  equal,  but  nonzero, 
magnitude  with  expectation  value  given  by  Eq.  (6.3).  When,  in  addition,  =  0  the 
signal  pattern  will  simply  be  the  complex  sum  of  the  true  and  image  fn  quency 
patterns,  as  in  Figure  14  for  the  quiescent  case,  with  a  null  in  the  jammer 
direction. 

Figure  17  shows  simulated  patterns  for  configuration  B  when  using  Applebaum’s 
algorithm  when  t Oj  =  0  and  Wy  *  0. 

For  the  Frost  case  with  a) i.  =  0  but  (O^  *  0  the  image  frequency  signal  -0)^  will  still 
act  as  a  new  jammer  giving  an  increased  sidelobe  level  coimibution  given  by 
Eq.  (6.5)  with  p  replaced  by  p,,  giving  -34  dB  for  configuration  B.  When  co  *  0  but 
ou^  =  0  the  desired  signal  will  be  deformed  as  in  Section  5,  resulting  in  a  very 
deformed  pattern.  The  situation  ov  =  toy  =  0  corresponds  to  a  complete  correlation 
between  the  desired  signal  and  jammer,  a  situation  which  will  always  give  very  bad 
patterns. 

Figure  18  shows  simulated  patterns  lor  configuration  B  wt-en  using  Frost's 
algorithm  when  oij  =  0,  (0^  *  0  and  when  w  *  0,  u)^  =  0. 
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Figure  17.  Adapted  Patterns  for  Configuration  B,  Applebaum’s  Algorithm. 
Uncorrelated  errors  a  =  p  =  0.07  (RMS).  (Dy=0,  co jk).  -  true  frequency, 

- image  frequency,  . -  no  IQ  errors. 


Relative  gain  [dB] 


6.4  Pattern  Effects  for  Identical  Channels 

In  Figure  19,  the  same  kind  of  patterns  as  in  Figures  15  and  17  for  Applebaum’s 
algorithm  are  shown  when  the  errors  are  identical  in  all  channels,  a  =  P  =  0.07, 
giving  a  sidelobe  peak  of  -26  dB  from  Eq.  (6.4).  Due  to  the  change  in  the  sign  of 
the  phase  of  the  error  signal,  a  signal  from  direction  0  will  generate  error  signals 
that  will  look  like  a  plane  wave  from  direction  -0.  To  illustrate  this,  the  main 
lobe  has  been  moved  to  10°  and  the  jammer  to  31°.  The  image  frequency  peak  at  0  = 

O  c 

-Qj  =  -10  is  given  by  Eq.  (6.4).  There  is  also  an  extra  null  at  -31  .  This  happens 
because  the  patterns  for  the  two  frequencies  are  mirror  images  of  each  other. 

Therefore,  generating  a  null  at  31°  for  the  image  frequency  will  cause  a  null  for 
the  true  frequency  at  -31°,  and  vice  versa.  Two  degrees  of  freedom  are  still 
required. 


7.  DC  OFFSET  ERRORS 


The  DC  offset  errors  in  the  video  circuits  and  A/D  converters  will,  as  opposed 
to  the  other  error  types  considered,  give  error  signals  that  are  independent  of  the 
input  signal.  Like  the  gain,  phase  shift,  and  IQ  errors  in  Sections  5  and  6,  it  may 
be  corrected  for  by  a  calibration  routine.  For  non-adaptive  antennas,  random 
DC-offset  errors  will  give  an  RMS  output  signal  level,  at  video  frequency  to  =  0, 
relative  to  the  maximum,  or  saturated,  power  of: 


DCoff 


<AZ> 

-  8  xr4-2-2b 

7nr 


(7.1) 


where  A  is  the  ratio  of  the  DC  offsets  (A)  and  the  quantization  step  ( q )  in  the  A/D 
converters,  <A^>  is  the  mean  square  and  b  is  the  number  of  bits  in  the  A/D 
converter.  Alternatively  we  can  write  the  DC  offset  signal  to  quantization  noise 
ratio  as: 


=  12<A2> 


(7.2) 


quant 
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Figure  19.  Adapted  Patterns  for  Configuration  A,  Applebaum’s  Algorithm.  Identical 
errors  a  =  (3  =  0.07. -  total  pattern  (u) j=0), - true  frequency, 

- .  image  frequency. 
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If  the  errors  are  correlated  between  the  channels,  or  even  identical,  the  DC 
offset  will  look  like  an  incident  plane  wave,  at  broadside  for  identical  errors, 
with  a  video  signal  frequency  of  zero  (center  of  RF  band).  However,  when  there  are 
large  gain,  phase  shift,  and/or  IQ  errors  that  are  compensated  for,  this 
compensation  will  tend  to  decorrelate  the  DC  offset  signals. 

In  adaptive  beamforming  the  DC-offset  signal  will  be  considered  as  a  jammer 
signal  by  the  process  and  will  thus  be  partly  eliminated.  This  will  require  one 
degree  of  freedom. 

Figure  20  shows  adapted  patterns  for  configuration  B,  using  Applebaum’s 
criterion,  when  the  RMS  value  of  the  offset  is  equal  to  the  RMS  noise  voltage  at  the 
input  of  the  A/D  converters,  giving  an  equivalent  total  jammer  power  of  P.  =  A/a2. 

It  is  here  assumed  that  the  video  frequency  (0  *  0  so  that  the  DC  offset  is 
uncorrelated  with  all  signals.  Patterns  for  both  uncorrelated  and  identical  errors 
are  shown  as  well  as  the  patterns  without  any  error.  Here,  the  desired  signal  is  at 

o  o  o 

20  and  the  jammer  at  43  to  better  show  the  null  at  0  for  the  pattern  with 
identical  errors.  This  null  may  give  a  severe  effect  when  the  desired  signal  is 
close  to  broadside.  It  is  seen  from  Eq.  (2.13)  with  P j/o2  =  N-A2/o2  that,  for  large 
N,  A  must  be  very  small  for  the  corresponding  null  depth  to  be  small. 

For  uncorrelated  errors  there  is  no  null  but  an  increase  in  sidelobe  level. 

Using  Eq.(2.I7)  with  Pj/o2  =  N<A2>/g2,  where  <A2>  is  the  mean  square  of  the  DC 
offset  errors,  we  get  the  average  sidelobe  level  contribution: 


Sidelobe  level  = 


[n  ♦  -^-1 

L  <A2>J 


(7.3) 


giving  -23.4  dB  in  the  simulated  case.  This  has  a  standard  deviation  equal  to  the 
average  itself  and  the  sidelobes  can  thus  vary  considerably  between  different 
simulations.  In  this  -40  dB  sidelobe  configuration  we  would  thus  need  at  least 
A  <  0.1a  when  the  errors  are  uncorrelated. 

RMS  v 

Since  no  additional  error  signals  are  introduced  by  the  desired  signal  the 
result  for  the  Frost  algorithm  is  the  same. 

For  digital  IQ-generation  there  is  no  DC-offset  error  for  the  Q  signal  if  this 
is  created  in  a  Hilbert  filter,  slopping  any  DC  component.  There  may  be  an  offset 
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error  for  the  I  signal  (possibly  converted  to  the  frequency  -/ JA)  but  as  there 
should  never  be  any  DC  component  in  the  signals  this  may  be  filtered  out  if 
necessary. 


-  No  errors,  —  uncorrelated,  identical 


Figure  20.  Configuration  B  with  DC-offset  Errors,  Using  Applebaum’s  Algorithm. 
- A=0, . uncorrelated  with  A  =o  ,  -  identical  with  A=a  . 

RMS  v  v 
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8.  SAMPLED  COVARIANCE  MATRIX 


In  the  previous  sections  it  was  assumed  that  the  covariance  matrix  used  was  the 
correct  one.  In  practice  we  don’t  know  the  expectation  values  of  the  matrix 
elements,  E[x.x*j ),  but  only  their  estimation  in  the  form  of  mean  values  over  a 
finite  number  of  time  samples,  Eq.  (2.8).  The  finite  number  of  samples  also  means 
that  some  nice  properties  of  the  covariance  matrix  used  in  Section  2,  such  as  no 
correlation  between  jammer  signals  and  noise  and  between  noise  in  different 
channels,  are  only  approximately  true. 


8.1  SNIR  as  a  Function  of  the  Number  of  Samples 


In  Reed15  it  is  shown  (see  also  Reference  3,  Chapter  6)  that,  for 
Applebaum’s  criterion,  K  >  2N-3  samples  are  necessary  for  the  expected  SNIR  to  be 
within  3  dB  of  its  asymptotic  value.  For  Frost’s  criterion,  more  samples,  or 
K  >  2A[1  +  P y(2o2)] ,  are  required  (Reference  3,  Chapter  6). 

We  can  define  three  different  types  of  SNIR’s  when  using  sampled  covariance 
matrices,  d>.  For  Applebaum’s  algorithm  the  "maximized"  SNIR  is  the  value  of  the 
maximized  function,  Eq.  (2.1).  When  x s^: 


SNIR  = 

max 


T  +  * 

w  xdxd  w 
wTd>  w* 


. .  I* 

sri>  s, 
a  a  a 


(8.1) 


which  is  equal  to  Applebaum’s  [4,  Eq.  (2.43)]  in  the  asymptotic  case.  This  is 
approximately  related  to  the  asymptotic  value  SNIR^m,  obtained  with  an  infinite 
number  of  samples  (using  the  correct  asymptotic  covariance  matrix)  by: 


SNIR  »  SNIR 

max  asym 


l 


FT 

~T 


(8.2) 


where  K  is  the  number  of  time  samples  used  in  the  estimation  of  <J>a< 


15  Reed,  I.S.,  Mallett,  J.D.,  and  Brennan,  L.E.,  (1974),  Rapid  convergence  in 
adaptive  arrays,  IEEE  Trans.  Aerospace  and  Electronic  Systems ,  Nov.,  pp.  853-863. 
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The  reason  that  SNIR  >  SNIR  is  that  the  sampled  noise,  together  with  the 

max  asym 

jammers,  can  be  considered  as  K  jammers  (one  strong  and  K-l  weak  when  there  is  one 
real  jammer).  Of  these  jammers  N- 1  can  be  canceled  leaving  K-(N-l)  weak 
noise-jammers.  [Eq.  (8.2)  can  be  shown  approximately  by  an  series  expansion, 

Eq.  (8.8),  of  <J>'‘  and  it  has  the  correct  behavior  when  K  =  AM  and  when  N  =  1.] 

More  important  is  the  expected  value  of  the  SNIR,  obtained  when  using  the 

A  A 

calculated  weights  and  random  uncorrelated  noise.  This  is.  when  x s^,  given 
by,  (Reference  15,  Eqs.  (15)  and  (18)): 


SNIR  = 

exp 


w" xdxd  w‘ 
wTd>  w* 

A 


SNIR  fl- 

asym  I 


N- 1  j 

rarj 


(8.3) 


8.2  No  Signal  Incident  on  the  Array 


This  is  the  situation  for  the  Applebaum  algorithm  with  no  jammer  and  for  the 
Frost  algorithm  with  no  jammer  and  a  very  weak  desired  signal.  The  effects  of 
sampling  will  enter  as  errors  in  the  estimation  of  the  noise  covariance  matrix.  This 
estimation  error  is  the  essential  factor  determining  the  resulting  sidelobe  level 
and  gain  loss  in  the  Applebaum  case  and  the  absolute  sidelobe  level  in  the  Frost 
case.  This  is  also  true  when  there  are  jammers  present  except,  of  course,  close  to 
the  jammer  directions. 

That  K  is  large  enough  to  give  a  good  expected  SNIR  as  discussed  above  does  not 
necessarily  mean  that  the  pattern  looks  good  in  terms  of  average  sidelobe  levels, 
etc. 

If  we  estimate  the  elements  of  the  noise  covariance  matrix  using  K  samples: 


<t> 


mj 


=  > 
‘  J 


i  ^ 

K  l 

k=l 


x-  ,x-  , 
i.k  j,k 


(8.4) 


the  following  statistical  relations  can  be  shown,  assuming  that  x.  ^ j  and  j c .  ^  are 
uncorrelated  when  i  *  j  or  kl  *  k2: 

Expectation  value:  E(d>^. )=  o28-j  (8.5) 
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so  that: 


Standard  deviation: 


/  Varj  l$n  J  ) 


VJC 


(8.6) 


$  =  o2ri  +  ir 

n  (  yr  J 


(8.7) 


where  a2  is  the  noise  power  in  each  channel,  I  is  the  identity  matrix  and  R  is  a 
Hermitian  matrix  where  each  element  is  a  random  variable  with  expectation  value  =  0 
and  standard  deviation  =  1,  (equally  distributed  in  the  real  and  imaginary  parts  for 
the  off-diagonal  elements).  As  an  example,  Figure  21  shows  results  from  a 
simulation  with  two  complex,  normally  distributed,  random  variables  with  cr2  =1.  The 
upper  left  figure  shows  the  magnitude  of  the  mean  values  as  function  of  the  number 
of  samples,  whereas  the  upper  right  shows  the  magnitude  of  the  variances  or  the  mean 
square  sum,  and  ^>m2<  and  the  lower  left  shows  the  magnitude  of  their 
covariance,  14^  I.  Lines  for  the  expected  values  and  the  standard  deviations  are 
also  shown. 

As  seen,  the  elements  in  the  covariance  matrix  approach  their  asymptotic  values 
rather  slowly.  If  VK~'  is  large  enough  so  that  the  eigenvalues  of  R/VJC  ,  which  are 
of  the  order  VNJK,  are  less  than  1,  the  inverse  can  be  found  using  the  series 
expansion  (Reference  2,  Appendix  1): 


[I  +  A]'1  =1  -  A  +  A2  -  A3  +  A4-.... 


(8.8) 


This  gives  the  inverse  of  the  covariance  matrix  to  the  first  order: 


(8.9) 


and  the  resulting  complex  conjugate  of  the  quiescent  Applebaum  weight  vector: 


•  •  i  *  i  _ 

w  =  w - R  w  -  s,  -  —  Rs  , 

q  qa  VK-'  qa  a  /IT  a 


(8.10) 


where  w  is  the  asymptotic  quiescent  weight  vector.  Thus  the  weight  vector  w 
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Covariance  Mean  value 


Samples 


Figure  21.  Two,  Normally  Distributed,  Simulated  Sequences,  a)  Mean  values, 
b)  Variances,  c)  Covariance.  -  -  -  -  Expectation  value  and  standard  deviations. 
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will,  to  the  first  order,  have  a  normalized  mean  square  error,  Eq,  (A.3)  of: 


<A2>  =  N/K  (8.11) 

which  will  result  in  a  mean  sidelobe  level,  Eq.  (A.6),  of: 

<A2>  , 

Sidelobe  level  =  ~j^r =  ‘  (8.12) 

'a  'a 

The  effects  of  the  random  matrix  elements  can  be  decreased  by  adding  an 
artificial  noise  covariance  matrix  c  I,  Eq.  (2.40),  increasing  the  magnitude  of  the 

art 

diagonal  elements.  This  will  make  the  relative  standard  deviation  of  the  matrix 
elements  a  factor  (1  +  o^o2)  smaller.  The  results  for  the  normalized  errors  and 
sidelobe  levels  will  still  apply  if  we  simply  replace  K  by: 

K  =  K(  1  +  a2  /a2)2  (8.13) 

art  arr 

Since  the  SNIR  when  there  are  no  jammers  is  equal  to  the  signal  power  per 
channel  times  the  antenna  gain  we  can  combine  the  rigorous  result  in  Eq.  (8.3)  with 
the  gain  loss  in  Eq.  (A.2)  to  get  a  higher  order,  "rigorous",  weight  error.  (The 
result  in  Eq.  (8.3)  was  derived  in  Reference  15  for  or  rlA(0j)  =  1.)  The 

result  of  this  is: 


<A2> 

R 


N 

K  +~\~~N 


<A2> 

_ i 

1  -  <A2> 
i 


(8.14) 


The  antenna  aperture  efficiencies  obtained  from  several  simulations  are  shown  in 
Figure  30  together  with  the  efficiency  obtained  from  Eqs.  (A.2)  and  (8.14): 

r  N  -  l/rj0(9  .)  , 

G(Qd )  =  N  T}(Qd)  =  Mye^l - K  T  ;  ■  -  j  (8.15) 

where  G(Qd)  is  the  antenna  gain,  *1(6^)  is  the  resulting  aperture  efficiency  and 
V9  d)  is  the  asymptotic  efficiency. 
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The  absolute,  sampling  induced,  sidelobe  gain  is,  from  Eq.  (A.2),  (see  also 
Reference  8): 


N 


Sidelobe  gain  =  -y-  — 


and  the  relative  sidelobe  level  is,  from  Eq.  (A.5): 


(8.16) 


Sidelobe  level  = 


r\{bdH«+l)  ~  vQLBjflK+l-N) 


TT 


(8.17) 


This  and  the  results  from  several  simulations  are  shown  in  Figure  31. 

Figure  22  shows  adapted  patterns  for  basic  configuration  B  using  Applebaum’s 
algorithm  and  K  =  32  (=2 N),  128,  512  and  2048  samples  together  with  the  asymptotic 
{K  =  °°)  result.  As  seen,  the  mean  sidelobe  level  is  high  for  the  lower  K.  Equation 
(8.17)  gives  sidelobe  contributions  of  -11,  -19,  -26,  and  -32  dB  respectively. 

Figure  23  shows  the  results  with  K= 32  when  using  o2  =  0,  10  a2  and  100  o2 

an 

respectively  together  with  the  asymptotic  pattern.  The  sidelobe  contributions  from 
Eq.  (8.17)  with  the  modification  of  Eq.  (8.13)  are  -11,  -35,  and  -54  dB  respectively. 

The  Frost  case  with  a  very  weak  signal  is,  as  was  shown  in  Section  4,  very 
similar  to  the  Applebaum  case  with  the  difference  that  the  tapered  quiescent 
weights  can  be  obtained  from  a  modification  of  the  covariance  matrix,  Eq.  (4.2), 
rather  than  from  a  tapered  steering  vector. 

Applying  Eq.  (4.2)  to  Eq.  (8.7)  we  get  to  the  first  order: 


= 

/IF 


D  1  [« _ o2  D  1  R1 

t2  L  a2  v'tH 

art  a  rt 


(8.18) 


to  be  compared  with  Eq.  (8.9)  for  the  Applebaum  case.  The  gain  and  sidelobe  level 
results  for  the  Applebaum  case  above  can  then  still  be  used  if  K  is  replaced  by: 


<w2>o4 

D 


(8.19) 


where  <w2>  is  the  mean  square  of  the  elements  in  D1. 
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Angle  fdeg] 


Figure  23.  Configuration  B.  Applebaum's  Algorithm.  Asymptotic  (A=°o)  and  K 
with  added  artificial  noise  cr  =  0,  10  a2  and  100  a2 


The  reason  for  the  pattern  degradation  in  the  basic  algorithm  is  the  random 
matrix  R,  which  makes  the  noise  eigenvalues  nonequal.  One  method  to  partially  solve 
the  problem  is  to  calculate  the  eigenvalues  and  eigenvectors  and  try  to  identify  the 
noise  eigenvalues,  which  are  of  the  order  of  o2  with  standard  deviation  a2VNJK,  and 
set  them  to  their  "correct"  value  o  ,  or  to  some  other  constant  value,  Eq.  (2.39). 

As  the  noise  eigenvectors  only  have  to  span  the  noise  space  they  can  still  be  used 
as  they  are. 

Since  the  signal  eigenvalues  are  often  significantly  larger  than  the  noise 
eigenvalues  the  identification  of  the  eigenvalues  is  often  not  very  difficult.  When 
the  jammers  are  so  weak  that  this  is  not  the  case  they  are  probably  not  so  important 
to  cancel  out.  With  this  method,  the  error  induced  sidelobes  will  be  very  small  and 
the  gain  is  essentially  equal  to  the  error  free  gain. 

Figure  24  shows  adapted  patterns  using  K-  32  samples,  when  setting  all 
eigenvalues  less  than  3a  equal  to  a  ,  (the  largest  noise  eigenvalue  in  this  example 
was  2.1a  ).  As  seen  there  is  a  big  improvement  in  the  sidelobe  region.  This  method 
is  also  discussed  in  Reference  7. 

8.3  The  Applebaum  Case  with  One  Jammer 


Now  the  sampled  covariance  can  be  written: 


<f>  =  x  x+ .  +  <I>  .  =  (x.  +  n„)(x.  +  n„)+  +  n,n! 

A  IS  t  IS  i  M  N  f  tl/v  i  7  1  l 


(8.20) 


where  x  .  is  a  "total '  jammer  vector,  being  the  sum  of  the  true  jammer  vector,  x ., 

K7  J 

and  a  noise  vector,  n|( ,  containing  the  time  correlation  between  the  K  noise  samples 

in  each  channel  and  the  jammer.  The  noise  vector,  nx,  is  the  rest  of  the  true  noise 

vector,  x^=  n(|  +  nx,  and  is  uncorrelated  with  \..  The  matrix  n±n^  is  equal  to 

the  sampled  noise  covariance  matrix,  <t>^,  to  the  order  o2//!^.  (The  elements  are 

slightly  smaller  however,  giving  a  slightly  smaller  gain  decrease  so  that,  including 

the  jammer  signal,  the  expected  SNIR  is  still  given  by  Eq.  (8.3).)  The  elements  of  n|( 

will  have  expectation  value  zero  and  standard  deviation  a/VK^.  Using  Eq.  (2.10)  with 

A  =  or  Eq.  (2.27),  the  weights  for  one  jammer  are: 
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Relative  gain  [dB] 


Modification  of  eigenvalues 


Angle  [deg] 


Figure  24.  Patterns  for  Configuration  B,  Using  Applebaum’s  Algorithm.  K  =  32. 
-  All  eigenvalues  less  than  3o2  set  to  o'.  -  -  -  Asymptotic  pattern. 


(8.21) 


2 

where  P  =  P  ■  +  Na  IK  is  the  total  power  in  the  total  jammer  vector  x  ..  This  gives  a 

k  j  KJ 

quiescent  pattern,  with  sampling  induced  sidelobes  and  a  "jammer  plus  correlated 
noise"  cancellation  beam. 

If  we  assume  that  there  is  a  perfect  null  in  the  x  .  =  x.  +  n„  "direction",  the 
null  depth  in  the  true  jammer  direction,  Xy,  is  obtained  by  projecting  the  n|(  vector 
onto  the  w  vector.  Assuming  n(|  is  a  random  vector  in  element  space,  the  expected 
result  is: 


Jammer  gain 


N  a1  a2 
-KIT-  Tp-j 


(8.22) 


where  Pj  is  the  jammer  power  per  channel.  This  jammer  gain  has  a  standard  deviation 
approximately  equal  to  the  expectation  value  itself.  The  null  depth  is  obtained 
from  Eqs.  (8.22)  and  (8.15). 


G(0.) 

Null  depth  =  ^  = 


(8.23) 


This  is  illustrated  in  Figure  32. 

For  more  than  one  jammer  the  result  is  essentially  the  same  because  nulling  any 
linear  transformation  of  the  jammers  (to  get  uncorrelated  jammers)  will  also  null 
each  separate  jammer. 

In  the  simulations  shown  in  Figure  22  the  null  depth  for  K  =  32  =  2N  was  -67  dB 
while  Eq.  (8.23)  gives  -63  dB. 

When  artificial  noise  is  added  Eqs.  (8.22)  and  (8.23)  are  still  valid,  using  the 
true  a2  and  K,  but  the  aperture  efficiency,  T](0^),  is  now  obtained  from  Eq.  (8.15) 
using  K  from  Eq.  (8.13).  In  the  simulations  in  Figure  23  the  null  depths  are  -73, 

-71,  and  -70  dB  respectively  while  Eq.  (8.23)  gives  -66  to  -63  dB.  This  is  also 
illustrated  in  Figure  22. 

In  the  patterns  in  Figure  24  where  the  noise  eigenvalue  modification  has  been 
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used,  the  absolute  jammer  gain  is  essentially  the  same  as  in  Figure  22.  However, 
since  the  array  gain  has  increased,  the  null  depth,  Eq.  (8.23)  has  improved.  This  is 
further  illustrated  in  Figure  32. 

When  using  this  eigenvalue  modification  method,  even  fewer  than  K  =  2N  samples 
may  be  used  because  the  random  matrix  R  is,  to  a  large  extent,  eliminated  and  the 
sampling  only  enters  through  the  projected  noise  vector  n((.  Figure  25  shows 
patterns  for  configuration  B  using  Applebaum’s  algorithm,  with  K  =  32,  16,  and  8 
samples  together  with  the  asymptotic,  K  =  pattern.  As  seen,  all  patterns  have 
good  sidelobes.  The  null  depths  have,  however,  decreased  to  -71,  -62,  and  -55  dB  in 
these  simulations  as  compared  to  the  asymptotic  -144  dB  and  -66,  -63,  and  -60  dB 
from  Eq.  (8.23)  with  ^(0^)  being  the  error  free  efficiency. 

In  the  basic  algorithm  the  reason  for  the  loss  of  expected  SNIR  is  mainly  loss 
of  antenna  gain,  while  the  null  depths  are  usually  good.  When  using  the  noise 
eigenvalue  modification,  and  when  adding  artificial  noise  with  sufficient  magnitude, 
the  antenna  gain  is  very  close  to  the  error  free  gain  and  the  reason  for  the  loss  of 
SNIR  is  the  decreased  null  depth,  increasing  the  received  jammer  signal.  Using 
Eq.  (8.22)  for  the  expected  jammer  gain,  the  expected  SNIR  becomes  for  these  methods, 
when  there  is  one  jammer  present,  approximately: 

SNIR  (G=G  )  «  SNIR  .  (8.24) 

exp  asym  asym  1  +  A 

In  the  simulations  in  Figure  25  the  expected  SNIR  decreases  to  20.8,  20.5,  and 

19.4  dB  compared  to  the  asymptotic  20.9  dB  and  20.7,  20.6,  and  20.3  dB  from 

Eq.  (8.24).  Like  the  null  depth,  the  SNIR  has  a  large  standard  deviation  so  that  for 
few  samples  these  numbers  may  vary  significantly  from  simulation  to  simulation. 

Simulated  SNIR  results  are  illustrated  in  Figure  33. 

8.4  The  Frost  Case  with  Desired  Signal  But  No  Jammer 

Here  the  situation  is  similar  to  the  Applebaum  case  with  one  jammer.  The  total 
sampled  signal  vector  x^  =  x^  +  n(|  is  now  nearly,  but  not  exactly,  parallel  to  the 
true  signal  vector  x^  (which  is  here  assumed  equal  to  the  steering  vector  s^)  so 
that  a  null  will  not  result  but  the  array  gain  may  be  decreased  significantly  for 
small  K. 
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Relative  gain  [dB] 


No.  of  samples:  8,  16,  32,  inf. 


Figure  25.  Patterns  for  Configuration  B,  Using  Applebaum’s  Algorithm.  K  =  32,  16, 
8  and  «>.  All  eigenvalues  less  than  3cr  set  to  cr. 


The  random  part  of  the  remaining  covariance  matrix  =  n±n^  will  also  give  a 
gain  decrease  as  discussed  in  Section  8.2  and,  as  seen  from  Eq.  (8.21),  the  total 
aperture  efficiency  is  essentially  the  product  of  the  quiescent  efficiency,  affected 
by  the  errors  in  4>  j_  and  obtained  from  Eq.  (8.15),  and  an  efficiency  affected  by  the 
error,  n,,,  in  x^. 
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ar  a'Vi 


(8.25) 


Using  the  fact  that  is  a  random  vector  with  standard  deviation  a/'/FT  in  all 
elements  we  can  write  the  unit  total  signal  vector  *Kd=  x^/ 1  *Kd  I  to  second  order  in 
G/vTTFJ 


K  d 


=  X 


Cl  N-\  o2l 


+  r 


(8.26) 


where  r  is  a  random  vector  where  each  element,  except  along  x^,  has  expectation 
value  =  0  and  standard  deviation  =  1.  Since  the  errors  are  included  in  the  adaptive 
process  the  resulting  weight  errors  will  be  correlated  and  Eq.  (A.2)  cannot  be  used 
in  the  main  beam  region.  Instead  we  can,  if  we  consider  the  total  signal  x  ^  as  a 
jammer,  get  the  array  gain  in  the  true  signal  direction,  x^,  by  using  Eq.  (2.15). 

The  quiescent  aperture  efficiency,  t]  (0.),  in  the  "jammer"  direction,  x  ,,  is  now 
for  uncorrelated  channels  obtained  using  Eq.  (8.26).  When  x^  =  s^: 


V  =|vv/= 1  • 


TP/-" 


(8.27) 


so  that  from  Eq.  (2.15): 


Gf(fid)  =  Nx](Qd)  =  H\(Qd) 


0  + 


N-  1  ,2 

~i r  J 


1  +  nr(7^ 


(8.28) 


where  T]  (0  ,)  is  the  quiescent  efficiency  in  the  true  signal  direction,  dependent  on 
the  errors  in  <J>^  and  is  obtained  from  Eq.  (8.15).  Since  we  have  neglected  the 
effects  of  the  errors  in  <J>^X  on  rj(| ,  Eq.  (8.25),  this  should  be  considered  as  a  first 
order  approximation  for  Gp(Qd),  valid  for  not  too  small  K.  Also,  since  we  have 
approximated  the  expectation  value  of  a  ratio  by  the  ratio  of  the  expectation 
values,  we  can  expect  some  errors,  especially  for  small  K. 

When  the  covariance  matrix  modification  for  the  Frost  algorithm,  Eq.  (4.2),  is 
used,  we  get,  using  Eqs.  (4.3)  and  (8.26)  and  the  same  method  as  when  deriving 
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Eqs.  (2.15)  and  (8.28): 


G?(Qd)  =  Nr)(Qd)=  Nr)(Qd) 


(1  +  <wd>-^1  )2 


1  +  <w  > 

D 


1(<WD  >Pd/°2+2) 


(8.29) 


where  <wq>  is  the  mean  value  of  the  elements  in  D’1  and  where  T|  (0^)  is  obtained 
from  Eqs.  (8.15)  and  (8.19).  When  a2,  <wq>  should  be  replaced  by  <wd>  o2 /O1  • 
Simulated  and  calculated  curves  for  the  aperture  efficiency  tj(O^)  are  shown  in 
Figure  30. 

The  error  in  the  unit  signal  vector  will  also  give  an  increased  sidelobe  level. 

In  the  sidelobe  region  the  errors  will  add  incoherently  and  Eq.  (A.2)  can  be  used.  As 
above,  we  can  separate  the  weight  errors  into  those  due  to  an  incorrectly  sampled 
noise  covariance  matrix  and  those  due  to  an  incorrect  steering  vector  being 
nonparallel  to  the  total  signal  vector  x From  Eqs.  (2.24)  and  (2.25)  the  total 
weight  error  can  be  written  as: 


<A2>  =  <A2>  +  <A2>  +  <A2xA2>  (8.30) 

n  s  n  s 

where  <Afi>  is  due  to  the  error  in  <t>nl  and  is  given  by  Eq.  (8.14)  together  with 
Eq.  (8.19),  and  where  <A2>  is  due  to  the  steering  vector  error.  Using  Eq.  (8.26)  in 
Eq.  (2.25)  we  get: 


<A2>  =  14- 
s  a2 


(8.31) 


where  K  should  be  replaced  by  K  from  Eq.  (8.19)  if  the  covariance  matrix 
modification,  Eq.  (4.2)  is  used. 

The  relative  sidelobe  level  is  obtained  by  using  Eq.  (A. 7)  together  with 


Eqs.  (8.28)  or  (8.29).  This  is  illustrated  in  Figure  31. 

Figure  26  shows  patterns  for  configuration  B  with  o2  =  a2,  10a2  and  100a2.  The 

art 

diagonal  elements  are  given  by  D  =l/(Chebyshev  weights)  with  w  =1.  The  sidelobe 

ii  Di^ax  j 

level  contributions,  calculated  as  described  above,  are  -17  dB  for  a  =  10a  and 

art 

-37  dB  for  c2=  100o2 
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8.5  The  Frost  Case  with  Desired  Signal  and  Jammer 


Due  to  the  sampling,  both  the  thermal  noises  and  the  desired  signal  will  be 
partly  correlated  with  the  jammer,  and  the  covariance  matrix  can  be  written  as: 


3>  =  x,  -xt .  +  x,  ,x!  ,+  , 

kj  kj  kd  kd  n- <- 


=  <xy  +  V  V  <xy  +  \#»+  n«/+  <V  +  "n/V  +  ni,/+  "X  <8-32> 


Frost:  added  noise  1,  10  and  100 


Figure  26.  Patterns  for  Configuration  B,  Using  Frost’s  Algorithm.  Asymptotic  (K=°°) 
and  K  =  32  with  added  artificial  noise  o2  =  a2,  10a2,  and  1002a. 

an 
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In  Eq.  (8.32),  x .  and  n..  •  are  the  parts  of  x  ,  and  x  being  correlated  with  x..  They 
a  .  .  J  i  i  _  a  i  \n  _  J 


where  is  the  rest  of  x 


have  expected  magnitudes  ix^  i  =v'/5  JK  '  and  I  n1( .  I  =oVN/K 

being  approximately  equal  to  x^,  and  where  n|(^  is  the  thermal  noise  which  is 

correlated  with  x^  but  not  with  \j,  being  approximately  equal  to  n(j  in  Eq.  (8.20). 

The  total  jammer  vector  x  .  is  now  uncorrelated  with  the  total  signal  vector,  x^, 

and  with  the  noise  vector  n^.  The  total  signal  vector  is  also  parallel  to  the 

true  signal  vector  x^  with  essentially  the  same  magnitude.  Assuming  that  the  total 

jammer  vector  x  .  has  been  completely  nullified,  the  null  depth  in  the  true  jammer 
*  KJ 

direction  x.  relative  to  the  gain  in  the  desired  signal  direction  is  obtained  by 
projecting  the  vector  (x^  +n|(y)/|  xKy|  onto  the  w  vector.  Assuming  that  n)( .  is  a 
random  vector  in  element  space  and  that  P .  »  P  ,,  the  expected  result  is: 

G(0.) 


d' 


'  »  ^ 

Null  depth  =  = 


(w7+,’J 


(8.33) 


where  r|(9^)  is  the  aperture  efficiency,  obtained  from  Eq.  (8.28)  or  (8.29).  The  null 
depth  will  have  a  standard  deviation  approximately  equal  to  the  expected  value 
itself.  The  null  depth  obtained  from  Eq.  (8.33)  and  from  several  simulations  is 
illustrated  in  Figure  32. 

In  the  simulations  in  Figure  26  the  null  depths  were  -55,  -49  and  -49  dB 
respectively  while  Eq.  (8.33)  gives  -45  dB. 

When  the  error  induced  sidelobes  have  been  reduced  to  a  small  value  by  using  the 
artificial  noise  method  or  steering  vector  projection  discussed  below,  the  expected 
SNIR  is,  as  for  Applebaum’s  algorithm,  decreased  from  its  asymptotic  value  mainly  by 
the  decreased  null  depths.  Using  Eq.  (8.33)  for  the  null  depth  we  get  for  one 
jammer: 


SNIR  =  SNIR  - - -  (8.34) 

exp  asym  +  SNIR^^ 

where  T)o  is  the  asymptotic  aperture  efficiency  and  r)  is  the  true  aperture  efficiency 
obtained  from  Eq.  (8.29).  This  is  illustrated  in  Figure  33  together  with  simulated 
results. 
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8.6  Effects  of  Steering  Vector  Projection 

Since  the  problem  with  the  high  sidelobes  in  the  Frost  algorithm  is  that,  due  to 
the  correlation  between  the  signal  and  noise  samples,  the  steering  vector  has  noise 
eigenvector  components,  the  steering  vector  projection,  Eq.  (2.41),  as  used  for  some 
other  error  types  may  be  used  here  too.  The  sidelobes  will  then  be  similar  to  those 
of  the  Applebaum  algorithm  using  a  steering  vector  with  small  errors.  The  steering 
vector  components  in  each  of  the  noise  eigenvectors  is  random  with  an  expected  RMS 
value  of  a/y/KP  j  .  The  total  mean  square  error  in  s^  after  the  projection  will  then 
be: 


<A  >  -  O  -gp - <r  -jrp - (8.35) 

a  da 

where  M  is  the  number  of  signal  eigenvectors.  The  resulting  sidelobe  level 
contribution  for  small  errors  is  then: 


_2 

Sidelobe  level  =  -g—p — ~  (8.36) 

r d  'a 

a 

or  a  IP^  times  that  of  the  basic  Applebaum  algorithm,  Eq.  (8.12).  The  gain  and 
aperture  efficiency  is  now  given  by  Eq.  (A. 2): 


G(Qd)  =  NT](Bd)=  N^(Qd) 


1  + 


TO 


1  + 


No2 

'XTi 


(8.37) 


while  the  null  depth  is  still  given  by  Eq.  (8.33)  and  the  SNIR  by  Eq.  (8.34). 

Figure  27  shows  adapted  patterns  for  Frost’s  algorithm,  configuration  A,  using 
K=  32  samples.  The  weights  have  been  calculated  from  a  direct  inverse  of  the  sampled 
covariance  matrix,  Eq.  (2.6),  by  using  a  projected  steering  vector,  Eq.  (2.41),  and  by 
using  the  asymptotic  covariance  matrix.  As  seen,  the  pattern  for  the  direct 
covariance  matrix  inversion  is  very  bad  while  the  pattern  obtained  when  using  the 
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projected  steering  vector  is  much  better.  The  sidelobe  level  using  Eq.  (8.36)  is 
-37  dB.  The  null  depth  is  only  at  -49  dB  in  this  example  while  Eq.  (8.33)  gives 
-45  dB. 

In  Figure  28  this  steering  vector  projection  has  been  combined  with  the 
modification  of  the  covariance  matrix.  Section  4,  to  get  the  Chebyshev  quiescent 
pattern  of  configuration  B.  Thus  first  the  noise  space  components  have  been 
eliminated  from  the  steering  vector,  Eq.  (2.41).  Then  the  noise  eigenvalues  have 
been  set  to  a2,  Eq.  (2.39),  as  in  Figure  24  for  the  Applebaum  case,  and  the  corrected 
covariance  matrix  has  been  obtained  using  Eq.  (2.31).  A  modified  covariance  matrix 
is  then  formed  by  modifying  the  diagonal  elements  as  in  Eq.  (4.2)  with  a2  =  a2. 

an 

Finally  the  weights  are  obtained  from  a  direct  inversion  of  this  matrix,  Eq.  (2.2). 

The  expected  sidelobe  level  from  Eq.  (8.36)  is  -36  dB  and  the  null  depth,  Eq.  (8.33), 
is  the  same  as  in  Figure  27,  -45  dB.  Also  shown  are  the  pattern  for  Applebaum’s 
algorithm,  as  in  Figure  24,  the  asymptotic  Frost  pattern  and  the  Frost  pattern 
obtained  when  only  modifying  the  covariance  matrix,  Eq.  (4.2). 

This  method  may  also  be  used  when  there  are  amplitude  and  phase  errors.  Figure 
29  shows  the  same  type  of  pattern  as  Figure  28  (except  the  direct  inversion  pattern) 

O 

when  there  are  0.3  dB  and  2  RMS  errors  present. 

8.7  Pattern  Parameters  Versus  the  Number  of  Samples 

Figures  30  to  33  show  simulated  and  calculated  results  for  the  signal  direction 
aperture  efficiency  T|(0^),  the  average  sidelobe  level,  the  null  depth,  and  the 
resulting  SNIR  respectively  for  basic  configuration  B.  All  simulated  curves  are 
averages  over  100  simulations.  The  solid  lines  are  the  simulaied  results  while  the 
dashed  are  the  calculated  curves  for  Applebaum’s  algorithm  and  the  dash-dotted  the 
calculated  for  Frost’s  algorithm.  The  figures  shows  results  when  using  direct 
inversion  of  the  covariance  matrix,  when  using  artificial  added  noise,  o  =10  C  , 

art 

and  when  using  noise  eigenvalue  modification,  including  steering  vector  projection 
in  the  Frost  case. 

The  curves  are  denoted  with  an  A  or  a  F  for  the  Applebaum  and  Frost  algorithms 
respectively.  The  direct  inversion  method  is  indicated  by  a  d,  n  indicates  the 
artificial  noise  method,  m  indicates  the  noise  eigenvalue  modification  in  the 
Applebaum  algorithm  and  pm  indicates  the  projected  steering  vector  and  noise 
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Figure  27.  Patterns  for  Configuration  A,  Using  Frost’s  Algorithm,  K=  32. 

- - direct  inverse  of  the  sampled  covariance  matrix,  -  projected 

steering  vector,  -  asymptotic  pattern. 
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loop=100,  artnoise  =  10,  -  sim,  —  calc  App,  calc  Fro 


Figure  30.  Simulated  and  Calculated  Aperture  Efficiency  for  Config  ration  B  as 

Functions  of  the  Number  of  Samples.  See  text.  -  simulations,  -  -  - 

Applebaum,  -  -  -  Frost 
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Figure  31.  Simulated  and  Calculated  Sidelobe  Level  for  Configuration  B  as 
Functions  of  the  Number  of  Samples.  See  text.  simulations, 

Applebaum,  .  Frost 
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Figure  32.  Simulated  and  Calculated  Null  Depth  tor  Configuration  B  as  Functions  of 

the  Number  of  Samples.  See  text.  -  simulations,  -  -  -  Applebaum,  . 

Frost 
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Figure  33.  Simulated  and  Calculated  SNIR  for  Configuration  B  as  Functions  of  the 
Number  of  Samples.  See  text.  -  simulations, - Applebaum,  —  Frost 
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eigenvalue  modification  in  the  Frost  algorithm.  (In  some  cases  it  may  be  difficult 
to  see  the  calculated  curves  since  they  are  so  close  to  the  simulated  curves.) 

Since  the  sampled  covariance  matrix  is  singular  for  K  <  /V  the  basic  Eq.  (2.2) 
cannot  be  used  in  this  case.  If  the  modification  of  the  Frost  covariance  matrix, 

Eq.  (4.2),  is  used,  this  matrix  will  not  be  singular.  In  the  Applebaum  direct 
inversion  curve  a  pseudoinverse  has  been  used  in  this  case.  This  is  equivalent  to 
using  a  partial  steering  vector  projection  method  where  the  steering  vector  is 
projected  on  the  "non-zero  eigenvalue  subspace".  This  is  probably  the  the  reason 
for  the  bump  in  the  efficiency  and  SNIR  curves  below  K=\6  samples. 

In  the  efficiency  curves.  Figure  30,  the  calculated  curves  are  obtained  using 
Eqs.  (8.15)  and  (8.13)  for  the  Applebaum  case  (except  for  the  eigenvalue 
modification  case  where  the  gain  is  very  close  to  the  asymptotic  value)  and 
Eq.  (8.29),  and  (8.37)  for  the  Frost  case.  The  agreement  is  very  good  except  for  the 
Frost,  direct  inversion  curve  where  the  calculated  is  slightly  lower  than  the 
simulated.  Here  the  relative  variance  of  the  denominator,  Eq.  (8.29)  is  rather  large 
so  that  taking  the  ratio  of  the  expectation  values  of  the  numerator  and  the 
denominator  will  give  some  errors.  In  fact,  taking  the  inverse  of  the  average  of 
the  inverse  efficiency,  l/<l/r|>,  gives  a  better  agreement  with  the  calculated 
result. 

In  the  sidelobe  curves.  Figure  31,  the  calculated  curves  are  obtained  from 
Eqs.  (8.13),  (8.14)  and  (A. 5)  for  the  Applebaum  case  (except  for  the  eigenvalue 
modification  case  where  the  sampling  induced  sidelobes  are  very  small)  and 
Eqs.  (8.30),  (8.14),  (8T9),  (8.31)  and  (A.7)  for  the  Frost  case  together  with  the 
calculated  efficiency  above.  The  simulated  curves  show  the  total  average  sidelobe 

O 

levels  of  the  pattern  outside  9  =  ±17.5  ,  including  the  regular  sidelobes  and  nulls. 

The  simulated  results  will  thus  go  asymptotically  to  approximately  -43  dB.  The 
calculated  results  on  the  other  hand  only  include  the  sampling  induced  effects. 

In  the  null  depth  curves  shown  in  Figure  32,  the  calculated  curves  are  obtained 
from  Eq.  (8.23)  for  the  Applebaum  case  and  Eq.  (8.33)  for  the  Frost  case  together  with 
the  calculated  efficiency  above.  The  simulated  and  calculated  curves  are  in  general 
in  very  good  agreement.  It  should  be  remembered  that  the  standard  deviation  in  each 
of  the  100  simulations  are  equal  to  the  expected  value  so  that  the  standard 
deviation  of  the  average  is  about  0.4  dB. 
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The  SN1R  results  are  shown  in  Figure  33.  Here  the  calculated  results  are 
obtained  from  Eqs.  (8.3),  (8.24),  (8.15)  and  (8.13)  for  Applebaum’s  case  and 
Eqs.  (8.34),  (8.29)  and  (8.37)  for  the  Frost  case.  Since  the  resulting  aperture 
efficiency  is  the  most  important  parameter  for  the  direct  inversion  cases  the 
agreement  is  also  similar  in  these  cases.  For  the  artificial  noise,  projected 
steering  vector,  and  modified  noise  eigenvector  cases,  the  aperture  efficiency  is 
high  and  the  SNIR  is  determined  mainly  by  the  null  depths.  The  agreement  between 
the  simulated  and  the  calculated  curves  for  the  Applebaum  algorithm,  using  these 
methods,  is  very  good  while  the  simulated  Frost  results  are  slightly  higher  than  the 
calculated  result  for  few  samples.  This  is  again  probably  mainly  an  effect  of  the 
approximation  of  taking  the  ratio  of  the  expectation  values.  As  seen  from  Figures 
30  and  32  both  the  gain  and  null  depths  are  well  approximated  with  these  methods. 

9.  SIGNAL  QUANTIZATION 


If  the  ratio  of  the  quantization  step  in  the  A/D  converters  (q)  to  noise  RMS 
voltage  (a  )  is  not  too  large  the  errors  due  to  signal  quantization  will  be  random 
and  uncorrelated  between  the  channels  (see  Figures  8  and  9  in  Reference  13).  This 
will  then  simply  increase  the  equivalent  noise  level  of  q  1 12  so  that, 
with  q  =  which  is  often  used,  the  results  with  q  -  0  are  applicable  if  we  put 
a2  =  —  a2  (0.35  dB  noise  increase). 

equiv  1 2 

Figure  34  shows  two  simulated  patterns  with  q  =  0  and  q  =  o  for  configuration 
A,  using  K  =  1024  samples,  which  are  very  similar.  In  these  simulations  the 
difference  in  the  resulting  SNIR  was  0.3  dB  and  the  null  depths  were  -88  and  -85  dB 
respectively  while  Eq.  (8.23)  gives  -82  dB  in  both  cases. 

The  small  quantization  steps  (one  q-step  in  one  1  or  Q  channel),  which  are 
present  in  the  total  signal  at  each  time  sample,  has  an  amplitude  of  V5/N  times  the 
quantization  noise  RMS  value  cr  V6/{  13.V)'  the  total  noise  RMS  amplitude  when  q  =  o  . 
When  many  samples  are  averaged  these  small  steps  will  disappear. 
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10.  WEIGHT  QUANTIZATION 


We  assume  here  that  the  weights  are  calculated  correctly  using  a  correct 
covariance  matrix  but  that  the  beamformer  can  only  use  the  first  6  bits  of  these 
weights.  There  is  one  sign  bit  and  the  least  significant  bit  =  2'(b'1).  These 
rounded  weights  will  not  really  be  the  optimal  6-bit  weights.  For  ordinary 
non-adaptive  arrays  the  corresponding  average  sidelobe  level  contribution  is: 

2  2  ’  2b 

Sidelobe  level  =  -» -  (10. 1) 

N<w2>r| 

where  <w  >  is  the  mean  square  of  the  weights  and  T)  is  the  aperture  efficiency.  As 
the  effects  of  the  rounding  were  not  considered  in  the  weight  calculation  we  can 
expect  the  effects  in  an  adaptive  array  to  be  similar,  including  the  null  depths. 

Like  the  finite  null  depth  due  to  sampling  in  Section  8,  the  null  depth  contribution 
in  Eq.  (10.1)  will  have  a  standard  deviation  equal  to  the  expectation  value  itself. 

Figure  35  shows  radiation  patterns  for  the  Chebyshev  standard  configuration  B, 
using  6,  8,  10  and  12  bits  plus  the  pattern  without  quantization.  It  is  seen  that 

O 

most  patterns  (except  for  6  bits)  are  very  similar  except  at  the  null  at  20  . 

However  in  this  example,  the  weights  with  6  bits  are  identical  with  the  6-bit 
quiescent  weights.  The  null  depth  for  the  8,  10,  and  12  bits  are  here  equal  to  -67, 

-71,  and  -74  dB,  to  be  compared  with  -144  dB  for  the  pattern  without  quantization, 
and  with  the  -57,  -69,  and  -81  dB  average  level  from  Eq.  (10.1). 

If  we  have  maximum  bad  luck  however,  all  weights  will  have  their  maximum 
quantization  error,  one  half  quantization  step,  adding  coherently  in  a  jammer 
direction.  The  maximum  sidelobe  contribution  due  to  this  is  equal  to  6N  times  the 
average  level  given  by  Eq.  (10.1),  or  approximately  20  dB  above  the  average  level  in 
this  example. 
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Figure  34.  Simulated  Patterns  for  Configuration  A,  A  =  1024.  Effect  of  signal 
quantization,  q  -  0  and  q  =  G^. 


81 


Angle  [deg] 


Figure  35.  Effects  of  Weight  Quantization  on  the  Patterns.  Configuration  B; 
number  of  bits  =  6,  8,  10,  12,  and  <».  (The  two  last  patterns  are  coincident 
in  the  plot  except  at  the  null.) 


11.  SATURATION 


For  strong  input  signals  the  receiving  system  becomes  nonlinear.  Here  only  the 
effects  of  this  nonlinearity  on  the  adaptive  process  will  be  considered.  Generally, 
a  nonlinearity  will  generate  new  frequency  components  (intermodulation  products)  of 
the  type  /  =  nf  +  mf ' ,  where  and  f2  are  two  input  frequency  components.  Usually 
the  components  with  n  -  1,  m  =  2  and  n  =  2,  m  =  1  are  the  most  important  since  they 
often  are  within  the  signal  bandwidth.  For  the  low  frequency  circuits,  where  the 
relative  bandwidth  is  large,  and  for  the  A/D  converters,  after  wh;ch  there  is  no  or 
little  filtering  but  a  possible  folding  into  the  signal  frequency  band,  the 
components  (harmonics)  with  n  or  m  -  0  may  also  be  important. 

If  the  nonlinearity  is  before  the  splitting  of  the  I  and  Q  components  the  output 
signal,  including  harmonics,  will  be: 

a=  a]  +  juQ=  a^xpOur+cp)  +  u3exp(3(ja)tt-<p))  +  a5exp(5(jo>r+q>))+...  (11-1) 

where  it  has  been  assumed  that  the  nonlinearities  are  symmetrical  so  that  only  the 
odd  frequency  components  are  generated.  All  a  are  nonlinear  functions  of  input 
power.  For  not  too  high  input  power  the  power  in  the  mo-component  will  be 
proportional  to  (input  power)^.  When  m  *  0  the  result  be  similar.  As  the 
relative  phases  in  the  channels  for  the  intermodulation  products  are  different  from 
those  of  the  input  frequencies,  these  signals  will  appear  to  come  from  different 
directions. 

If  the  nonlinearity  is  after  the  splitting  of  the  1  and  Q  components,  for 

example  as  a  limited  dynamic  range  of  the  I  and  Q  A/D  converters,  all  harmonics  in 

the  Q  channel  are  delayed  1/4  period  of  the  basic  frequency  component  so  that  the 

O 

extra  90  phase  shift  will  also  get  multiplied  with  n.  This  is  equivalent  to  a 

change  in  sign  for  the  frequency  components  n  =  3,  7,  11,  15  etc.,  giving  tht  output 

signal: 

a  =  a  +  j a  -  a  exp(jwr+<p)  +  a  exp(-3(jcot+(p))  +  a  exp(5(jcoH-<p))+...  (11.2) 

c  I  Q  1  3  5 
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If  the  I  and  Q  channels  are  not  identical,  both  the  positive  and  negative 
frequency  components  must  be  used,  as  in  Section  6  where  the  component  n  =  -1  was 
generated.  That  is  also  true  if  the  nonlinearity  is  nonsymmetric  so  that  even 
frequency  components  are  also  generated.  The  IQ  errors  in  Section  6  have  much  in 
common  with  the  saturation  effects,  the  main  difference  being  the  number  of 
harmonics  and  signal  power  dependence. 

If  the  nonlinearities  are  to  have  negligible  influence  on  the  adaptive  process 
the  total  power  in  the  new  frequencies  should  not  be  higher  than  the  noise  power. 

For  a  Gaussian  distributed  noise  jammer  this  means  that  the  average  jammer  power  per 
channel  should  not  be  larger  than  about  II  dB  below  the  saturation  level  for  10  bit 
converters  with  a  quantization  step  equal  to  the  noise  RMS  voltage  (q  =  a  ),  and 
about  15  dB  below  for  18  bit  converters,  provided  that  the  the  jammer  is  really 
Gaussian  distributed  far  out  in  the  tail  of  the  distribution,  which  could  be 
questioned.  For  CW  jammers  the  power  limit  is  more  distinct;  the  power  should  be 
below  the  saturation  level. 

As  an  example,  Figure  36  shows  simulated  radiation  patterns  for  three  different 
jammer  input  power  levels  for  basic  configuration  A.  The  covariance  matrix  has  been 
estimated  from  K  =  4096  samples,  corresponding  to  an  RMS  sidelobe  level  of  -36  dB 
from  Eq.(8.12). 

It  is  assumed  that  all  channels  are  identical.  Radiation  pattern  means  the 
pattern  that  the  calculated  weights  will  produce  in  a  linear  system.  This  is  also 
equal  to  the  radiation  pattern  of  the  basic  signal  component,  when  the  dynamic 
ranges  of  the  channels  are  identical.  For  identical  channel.:  the  other  induced 
frequency  components  will  have  similar  patterns,  shifted  in  angle.  The 
nonlinearities  are  assumed  to  be  lumped  in  the  two  12  bit  I  and  Q  A/D  converters. 

The  quantization  step  is  equal  to  noise  RMS  voltage  (q  =  o  )  so  that,  for  a  CW 

jammer,  saturation  begins  at  a  jammer  to  noise  ratio  per  channel  of  JNR  =  63  dB.  The 

saturation  is  assumed  to  be  of  a  hard  limit  type,  so  that  the  output  is  linear 
up  to  the  saturation  level  and  constant  above  this  limit.  Patterns  are  shown  for  a 
Gaussian  noise  jammer  with  a  jammer  to  noise  ratio  per  channel  of  JNR  =  53  dB,  where 

there  were  no  saturated  samples  at  alt,  for  JNR  =  58  dB,  where  the  nonlinearities 

O 

has  started  to  show  up,  and  for  JNR  =  70  dB.  When  the  jammer  is  at  0y,  18  here. 
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Figure  36.  Saturated  Patterns,  Configuration  A,  Applebaum’s  Algorithm,  K  =  4096 

samples:  -  JNR  =  70  dB, - JNR  =  58  dB,  JNR  =  55  dB.  The  + 

marked  numbers  denote  harmonic  number  from  Eq.(Il,3). 


the  equivalent  directions,  0^,  of  the  nth  harmonic  can,  for  half  wavelength  element 
spacing,  be  found  from  ( k  is  an  integer  to  make  0^  real): 


sin  0  =  n  sin  0.  +  2k  (11.3) 

n  j 

If  the  channels  have  large  errors  that  are  corrected  for,  these  corrections  will 
leave  the  harmonics  uncorrelated  between  the  channels.  These  harmonics  will  then 
not  give  any  deep  nulls  but  an  increased  average  sidelobe  level,  as  in  Eq.  (2.17). 

The  resulting  SNIR  will  however  be  approximately  equal  to  that  for  identical 
channels. 

The  directions  obtained  from  Eq.  (11.3)  are  indicated  in  the  figure  and  it  is 
seen  that  they  correspond  to  nulls  in  the  JNR  =  70  dB  pattern  for  quite  high  values 
of  the  harmonic  number  n.  This  shows  how  the  adaptive  process  will  try  to  eliminate 
all  relevant  signals  to  get  a  high  total  SNIR.  This  will  of  course  cost  some 
degrees  of  freedom,  which  may  be  important,  especially  for  an  antenna  with  few 
elements.  In  his  example  with  N  =  16  elements  and  one  jammer  the  antenna  would  be 
able  to  null  out  harmonics  up  to  n  =  29.  As  the  power  for  higher  n  may  still  be 
quite  high  compared  to  the  noise  (for  very  high  power  the  saturated  signal  will 
become  a  square  wave  and  the  power  in  the  harmonics  will  go  as  1/n  )  there  will  be 
some  residual  jammer  power  decreasing  the  SNIR.  Of  course,  the  algorithm  does  not 
eliminate  the  lowest  order  harmonics,  leaving  the  other  harmonics  unaffected,  but 
will  affect  all  harmonics  to  give  the  lowest  total  SNIR.  For  the  three  different 
jammer  powers  in  this  example  the  deviations  cf  the  resulting  SNIR  from  the 
asymptotic  nonsaturated  values  were  0.3,  9.5,  and  32.1  dB  respectively.  These 
numbers  are  for  one  particular  statistical  outcome  of  the  simulations;  the  numbers 
vary  somewhat  between  different  simulations,  in  particular  for  the  JNR  =  58  dB  case 
where  there  are  relatively  few  saturated  samples  The  SNIRs  here  are  the  maximized 
SNIR  obtained  from  Eq.  (8.1)  and  not  the  asymptotic  SNIR  or  the  expected  SNIR  using 
the  calculated  weights.  The  difference  for  K  =  4096  samples  is  very  small,  however. 

For  larger  antennas  more  harmonics  can  be  suppressed.  The  deviations  from  the 
asymptotic  nonsaturated  cases  in  simulations  with  N  =  16,  32,  64,  and  128  elements 
and  JNR  =  70  dB  were  equal  to  32.1,  22.3,  6.7  and  0.8  dB  respectively. 


86 


There  is  also  a  possibility  that,  as  for  n=  13  in  this  example,  a  jammer  outside 
the  main  beam  will  generate  a  frequency  component  with  an  equivalent  direction 
within  the  main  beam,  which  may  have  serious  effect  on  the  SNIR  and  main  beam  shape. 

When  the  desired  signal  is  present  during  the  adaptation,  as  in  Frost’s 
algorithm,  the  corresponding  intermodulation  products  will  also  be  considered  as  new 
uncorrelated  jammers  by  the  adaptive  process. 

One  possible  method  to  improve  the  pattern  and  the  SNIR  would  be  to  use  only  the 
nonsaturated  signal  vectors  when  forming  the  covariance  matrix.  This  assumes  of 
course  that  the  same  is  also  done  for  the  total  output  signal  from  the  antenna. 

Hence  much  of  the  information  in  the  desired  signal  may  be  lost.  Using  only 
nonsaturated  signal  vectors  in  the  JNR  =  70  dB  case  would  cause  only  about  20 
percent  of  the  samples  to  be  used.  Even  fewer  samples  would  be  used  for  a  saturated 
CW  jammer. 

For  less  hard  saturations,  without  any  knee  in  the  transfer  function,  the 
magnitudes  of  the  higher  order  harmonics  will  be  smaller.  These  less  hard 
saturations  are  more  realistic  for  amplifiers  and  mixers.  A  non-linear  transfer 
function  of  the  type: 

1/ 

i  n 

1  +  (V  /V  ] 

i  n7  max*7 

has  also  been  used  in  the  simulations  (p  =  °°  corresponds  to  the  hard  limit  case). 

When  p  =  6  the  compression  at  V  =  V  is  1  dB.  For  N  =  16  and  JNR  =  70  dB  the 
deviations  of  the  resulting  SNIR  from  the  asymptotic  nonsaturated  value  were  for  p  = 

6,  12,  and  24  equal  to  18.2,  25.7,  and  29.7  dB  respectively,  to  be  compared  with 
32.1  dB  for  the  hard  limit  case. 


Vp 


(11.4) 


12.  CONCLUSIONS 


The  effect  of  various  imperfections  on  the  pattern  of  adaptive  arrays  has  been 
investigated.  The  algorithms  considered  were  Applebaum’s  "Maximum  signal  to 
interference  plus  noise  ratio"  and  Frost’s  "Constrained  power  minimization".  In 
these  algorithms  there  are  no  explicit  constraints  or  requirements  on  the  sidelobe 
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levels  included  and  hence  some  types  of  array  imperfections  can  cause  severe 
degradation,  in  particular  in  Frost’s  algorithm  where  the  desired  signal  is  directly 
included  in  the  adaptive  process. 

Three  types  of  solutions  to  these  problems  have  been  discussed:  to  make  the 
errors  small  enough  or  to  compensate  for  them  when  they  are  known;  to  reduce  their 
influence  by  adding  artificial  noise;  and  to  use  some  a  priori  knowledge  about  the 
covariance  matrix  and  its  eigenvalue  spectrum.  In  many  cases,  where  the  error  is 
due  to  an  incorrectly  estimated  covariance  matrix  or  steering  vector,  the  methods 
discussed  will  also  improve  the  SNIR. 

It  may  also  be  desirable  to  modify  the  correct  covariance  matrix  slightly  to 
obtain  some  desirable  pattern  properties  at  a  very  low  cost  in  terms  of  resulting 
SNIR  as  discussed  in  Sections  3  and  4. 

One  type  of  error  occurs  when  the  steering  vector  used  is  not  optimally  chosen. 
This  is  the  case  for  the  amplitude  and  phase  errors  in  Section  5  and  for  the  IQ 
errors  in  Section  6,  when  the  signal  video  frequency  U)^  =  0.  The  Frost  algorithm  is 
most  affected  by  steering  vector  errors  because  some  of  the  desired  signal,  which 
does  not  match  the  steering  vector,  is  considered  as  a  coherent  jammer  by  the 
process.  These  errors  can  in  principle  be  measured  and  corrected  for. 

Alternatively,  they  can  be  partly  corrected  for  by  decreasing  the  importance  of  the 
desired  signal  by  adding  artificial  noise  or  by  using  a  priori  knowledge  about  the 
noise  eigenvalues  and  signal  subspace  eigenvectors. 

Another  type  of  imperfection  is  that  some  additional  jammer  signals  may  be 
created  internally  in  the  array.  This  is  the  case  for  the  IQ  errors  in  Section  6, 
when  the  video  frequency  co  *  0,  for  the  DC  offset  errors  in  Section  7,  and  for  the 
saturation  effects  in  Section  11.  The  first  two  could  be  compensated  for  if  they 
are  known,  while  the  last  one  should  be  kept  small  by  having  sufficiently  large 
dynamic  range  in  the  receiver  modules  and  A/D  converters,  or  by  rejecting  saturated 
samples. 

A  third  error  type  is  due  to  quantization  of  the  signals,  discussed  in  Section 
9,  and  quantization  of  the  calculated  weights,  discussed  in  Section  10.  The  first 
will  give  an  additional  equivalent  noise  contribution  while  the  second  will  limit 
the  obtainable  null  depth.  The  latter  may  be  very  serious  and  a  sufficient  number 
of  bits  is  necessary. 
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If  the  noise  covariance  matrix  is  not  an  identity  matrix,  as  is  usually  assumed, 
large  effects  on  the  pattern  may  result  as  discussed  in  Section  3.  Forcing  it  to  be 
equal,  or  close  to,  an  identity  matrix  gives  better  sidelobes  at  the  cost  of  a 
slightly  lower  SNIR. 

As  the  covariance  matrix  can  only  be  estimated  by  a  time  average  it  will  have 
some  errors.  This  problem  is  discussed  in  Section  8.  In  particular  there  will  be 
apparent  correlations  between  truly  uncorrelated  signals,  including  the  noise.  This 
type  of  error  can  also  be  partly  corrected  for  by  using  a  priori  knowledge  about  the 
eigenvalue  spectrum  or  by  minimizing  the  relative  errors  by  adding  artificial  noise. 
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Appendix 

Error  Induced  Gain  Loss  and  Sidelobes 


For  an  array  antenna  the  gain  is,  for  isotropic  element  patterns,  given  by: 


2. 

|£  W.  exp(y'kr.)| 

G(Q)  =  -l  =  1  'N - —  =  N- ti  (0) 

yu  I2 

ih  ■ 


(A.l) 


where  the  w  are  the  total  element  weights,  including  possible  gain  and  phase  shift 
errors,  k  is  the  propagation  vector  and  r-  is  the  element  position  vector.  When  the 
weights  have  random  errors,  the  total  weight  vector  w  can  be  written  as  a  sum  of  a 
vector  wq,  parallel  to  the  error  free  weight  vector,  and  a  vector  Aw,  uncorrelated 
with  the  error  free  weight  vector  so  that  w  =  wq+Aw  and  E|w*Aw)  =  0.  If  we  also 
assume  that  the  weight  errors,  Aw,  are  mutually  uncorrelated,  the  expected  gain  in 
direction  0  can  be  written  as: 
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(A.2) 


where  <A  >  is  the  normalized  mean  square  error: 


<A2>  =  1*2L 


w 


2  N,  ,, 

i  =  l  'i=l 


(A.3) 


and  Tj  (0)  is  the  error  free  aperture  efficiency  as  a  function  of  the  angle  0.  In 

Eq.  (A.2)  we  have  also  approximated  the  expectation  value  of  the  ratio  with  the  ratio 

of  the  expectation  values.  This  will  often  give  a  result  that  is  slightly  too  low. 

For  very  large  errors  we  can  see  from  Eq.  (A.2)  that  G(0)  -»  1. 

The  expected,  error  induced,  sidelobe  level  relative  to  the  error  free  gain  in 
the  desired  signal  direction  is  thus: 


Sidelobe  level  = 


<A2> _ 

"WC  +  <aV> 


(A.4) 


and  the  sidelobe  level  relative  to  the  actual  signal  gain: 


Sidelobe  level  = 


<A2> 

*n0(V+<A2> 


(A. 5) 


When  the  errors  are  small,  or  if  we  know  that  the  error  weights  Aw.  do  not  give 
any  contribution  in  the  main  beam  direction,  we  can  use: 


Sidelobe  level  = 


<A2> 


(A.6) 
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Sometimes,  such  as  when  the  weights  will  tend  to  partly  cancel  the  desired 
signal  in  the  Frost  algorithm,  the  components  in  Aw  are  mutually  correlated  and 
Eq.  (A.2)  is  valid  only  in  the  sidelobe  region.  The  relative  sidelobe  level  can 
still  be  found  from  Eq.  (A.4)  if  Ar)o(0^)  is  replaced  with  the  actual  gain,  G(0j). 

Sidelobe  level  =  - -  (A.7) 

G(Bd){  1  +  <A2>) 

If  the  antenna  elements  have  an  element  gain  factor,  g(8),  this  is  simply 
inserted  in  the  numerator  of  Eqs.  (A.  1)  and  (A.2)  and  in  both  the  numerator  and 
denominator  of  Eqs.  (A.4)  to  (A.7). 
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